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ABSTRACT 


This  report  is  concerned  first  with  the  computer  deter- 
mination of  optimum  parameters  of  a tapped-delay  line 
canonic  model  of  third  order  interference  generation, 
and  second  with  the  use  of  this  model  in  simulating  re- 
ceiver response  to  waveforms. 


Factors  such  as  complexity  and  realism  of  the  inter- 
ference environment  are  used  to  establish  the  necessary 
cases  and  boundary  values  to  be  approximated  by  the 
model.  The  convergence  of  the  tapped-delay-line  model 
to  a frequency-polynomial  model,  and  the  rate  of  this 
convergence,  are  used  to  establish  limits  on  group  delay 
and  tap  spacing.  ,x  An  algorithm  is  developed  embodying 
these  limits,  that  recursively  searches  for  the  best 
group  delay  and  tap  spacing,  and  explicitly  determines 
the  best  tap  coefficients.  Sampled  data  techniques  are 
developed  for  generation  of  AM/FM  or  noiselike  signals 
and  ir.terferers,  for  generation  of  interference  products 
using  the  model,  for  simulation  of  the  equivalent  linear 
filter  of  tue  cascade  model,  and  for  calculation  of  dis- 
tortion in  phase  demodulator  or  amplitude-demodulator 
outputs.  The  functioning  of  the  computer  programs  is  dem- 
onstrated with  printouts  from  the  tests. 
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SECTION  1 


INTRODUCTION  AZ©  SUMMARY 

This  technical  report  for  tne  period  8 Decenter  1973  to  7 March 
1974  discusses  the  oonputer  determination  of  optiram  parameters  of  a 
tapped-delay  line  canonic  model  of  third  order  interference  generation, 
and  second,  (tfrith  the  use  of  this  model  in  simulating  receiver  response  to 
waveforms. 

1 . 1 Introduction 

In  earlier  studies  by  SIGNATRON  extensive  effort  has  been 
devoted  to  the  determination  of  the  response  of  nonlinear  cir- 
cuits to  multiple  sinusoidal  inputs  as  characterized  by  the  non- 
linear transfer  function  of  the  circuit.  The  present  effort  is 
concerned  with  the  development  of  canonic  models  that  will  per- 
mit easier  determination  of  the  network  response  to  modulated 
input  signals. 

1.1.1  Specific  Program  Objectives 

The  development  of  canonic  models  falls  naturally  into  a 
sequence  of  steps  which  form  the  specific  objectives  of  this 
program: 

a.  The  determination  of  analytically  tractable  approxima- 
tions to  the  nonlinear  circuit  response  that  are  partic- 
cularly  useful  for  the  small-percentage-bandwidth  sig- 
nals of  interest  in  communications. 

b.  Determining  the  minimum  number  of  parameters  necessary 
to  characterize  these  approximations. 

c.  Determine  the  minimum  number  and  most  effective  set  of 
measurements  that  will  permit  the  model  parameters  to 
be  extracted,  both  for  entire  receivers,  and  for  com- 
ponent simplifiers  and  mixers. 

d.  Verification  of  the  feasibility  of  the  measurement 
procedure . 
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e.  Determine  necessary  computer  programs  to  calculate 
the  model  parameters  either  from  measured  data  or 
from  analytic  predictions  of  circuit  response. 

f.  Determine  necessary  computer  programs  for  prediction 
cf  response  to  specific  modulated  input  waveforms 
using  the  measured/calculated  model  parameters. 

1 . 2 Summary  of  this  Report. 

This  report  is  concerned  entirely  with  items  (e)  and  (f)  in 
the  list  of  objectives:  the  determination  and  generation  of  com- 
puter program©  to  calculate  optimum  model  coefficients  and  simulate 
the  effect  on  waveforms  cf  the  receiver  being  modeled. 

Section  2 defines  the  objectives  of  the  computer  programs 
and  establishes  appropriate  limits  on  model  complexity,  interference 
cases,  and  model  parameters. 

Section  3 deals  with  the  calculation  of  optimum  model  param- 
eters including  tap  spacing,  group  delay  and  optimum  tap  coef- 
ficients. The  rate  of  approach  to  a frequency-polynomial  approxi- 
mation is  used  to  define  search  limits  for  tap  spacirg  while  the 
limiting  behavior  provides  a guideline  to  setting  search  limits 
for  group  delay.  An  algorithm  is  then  developed  that  uses  a re- 
cursive search  for  optimum  group  delay  and  tap  spacing  and  an 
explicit  optimization  of  the  tap  coefficients. 

Section  4 derives  properties  of  second-order  sampled  data 
Butterworth  filters  used  in  signal  filtering  and  noise  generation. 

It  is  useful  to  have  available  an  analytically  well-understood 
nonlinear  circuit  for  testing  program  routines.  Section  5 describes 
such  a circuit  and  the  associated  formulas  for  calculation  of 
and  transfer  functions. 
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In  Section  6 we  discuss  the  generation  of  signal  and  interfer- 
ence waveforms  which  may  be  sinusoidally  amplitude  modulated, 
phase  modulated,  or  both,  or  may  be  noiselike  in  character.  These 
are  available  as  an  alternative  to  sampled  data  tapes  derived 
from  actual  signal  sources. 

Computer  programs  have  been  successfully  written  that  embody 
all  of  the  features  described  in  Sections  2 to  6.  In  Section  7 
we  discuss  the  results  of  an  end-to-end  test  of  these  programs. 

1.3  Contributors 

The  work  reported  on  here  was  performed  by  L.H.Vears,  J.N. 
Pierce,  N. Johnson,  H.Gish  and  S.H.Richman.  This  report  was 
prepared  by  Ms.  Vears  and  Mr.  Pierce. 

1.4  Acknowledgments 

We  are  indebted  for  program  guidance  and  technica1  suggestions 
to  Mr.  John  F.  Spina  of  RADC  and  Prof.D.  Weiner  of  Syracuse  Uni- 
versity. 
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SECTION  2 

COMPUTATIONAL  OBJECTIVES 


The  bulk  of  the  work  discussed  in  this  report  relates  to 
the  generation  of  computer  programs  related  to  the  canonic 
modelling  effort.  In  this  section  we  discuss  the  objectives 
of  this  software  effort. 

2.1  General  Objectives 


Our  Technical  Report  #2  demonstrated  that,  at  the  present 
time,  the  scope  of  clonic  modelling  should  be  limited  to  the 
modelling  of  third-order  nonlinear  transfer  functions.  We 
further  established  that  for  computational  purposes  only  a very 
few  models  were  practical.  For  the  purposes  of  this  effort  we 
have  chosen  the  most  useful  of  these,  the  generalized  tapped- 
delay  line  model,  which  involves  complex-exponential  approxima- 
tions to  the  third-order  nonlinear  transfer  functions. 

The  software  necessary  to  use  canonic  modelling  must  provide 
the  following  capabilities: 

a)  A program  to  calculate  the  parameters  of  the  tapped- 
delay  line  model  either  from  measured  values  of  the 
transfer  function  or  from  transfer  function  values 
calculated  by  programs  which  analyze  the  nonlinear 
circuits. 

b)  A program  that  accepts  arbitrary  signal  and  inter- 
ference inputs  consistent  with  the  model  bandwidth 
and  sampling  rates,  and  generates  the  complex  en- 
velope of  the  corresponding  third-order  interference. 

c)  Routxnes  to  generate  realistic  signa?  and  inter- 
ference inputs  to  be  used  with  the  program  in  (b) 
above . 

d)  Routines  to  display  the  effect  of  the  interference 
on  the  baseband  output  of  a receiver's  demodulator. 
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2 . 2 Spectral  Structure  of  Signals  and  Interference 


To  put  some  structure  on  the  computer  modelling,  we  can  visu- 
alize the  RF  spectrum  as  consisting  of  a large  number  of  equally 
spaced  channels  with  a separation  of  W Hz  between  the  center  fre- 
quencies of  adjacent  channels.  This  structure  is,  in  fact,  quite 
typical  of  military  spectral  allocations.  The  same  number  W will 
also  be  roughly  equal  to  the  typical  signal  bandwidth,  and  typical 
receiver  IF  bandwidth,  if  we  take  these  bandwidths  to  be  defined 
by  the  (-20  dB)  or  (-60  dB)  points  on  the  spectra,  for  example. 
Since  the  objective  here  is  the  modelling  of  nonlinear  effects, 
it  is  an  adequate  approximation  to  equate  the  -3  dB  bandwidths 
to  W as  long  as  the  software  routines  avoid  any  linear  adjacent 
channel  interference  effects.  The  basic  framework  will  then  be 
taken  as  a desired  signal  at  the  tuning  frequency  v,  and  potential 
interferers  at  v + W,  v + 2W,  etc.,  all  with  equal  bandwidth  W, 
which  is  also  to  be  taken  as  the  IF  bandwidth. 


Now  let  vi'v2'  V3  ke  t^le  carrier  frequencies  of  the  three 
interfering  signals.  Then,  as  was  pointed  out  in  TR  #2,  the 
bandwidth  of  the  third-order  interaction  is  3W  so  that  inter- 
ference to  the  desired  signal  can  occur  if 


+ v2  - v3  « v.  (2-la) 

or 

V1  + '‘2  “ v3  = v --  W*  (2-lb) 


The  most  general  type  of  computer  modelling  would  then 

admit 


b) 


c) 


an  arbitrary  signal  modulation  at 
frequency  v. 


the  carrier 


three  independent  interference  modulations  at 
carrier  frequencies  v]/v2'vV 

the  interaction  carrier  frequency  located  at  either 
v or  one  of  the  adjacent  carrier  frequencies  v + W. 
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If  we  take  into  account,  however,  the  relative  importance  or  re- 
lative probability  of  the  various  types  of  interference,  the 
scope  of  the  modelling  can  be  reduced  with  no  loss  of  utility. 

We  will  now  develop  these  specific  restrictions.  Before  pro- 
ceeding to  this  it  is  helpful  to  review  some  results  from  TR  #2 
and  TR  #3. 


2 . 3 Use  of  Equivalent  Receiver 

In  Section  2.1  of  TR#3  we  introduced  the  concept  of  an 
equivalent  receiver.  We  repeat  the  relevant  part  of  Fig.  2.5 
of  that  report  as  Fig.  2.1  here.  The  essence  of  the  equivalent 
receiver  is  to  replace  the  distributed  (H^-  H^)  structure  of 
the  actual  receiver  with  a single  linear  filter  (with  transfer 
funct-  on  H(f))  following  a parallel  combination  of  a unit- 
gain  amplifier  and  a third-order  transfer  function 

This  equivalent  receiver  structure  also  forms  an  excellent 
signal  flow  chart  for  computer  simulation  in  that  the  possibly 
complicated  tapped-delay  line  structure  for  the  equivalent  IF 
filter  H(f)  can  be  applied  to  the  one-dimensional  output  of 
the  third-order  filter  f^)  rather  than  having  its  ef- 

fects incorporated  in  the  three-dimensional  tapped-delay  line 
structure  which  synthesizes  the  nonlinear  response. 

2 . 4 Formulas  for  Third-Order  Zonal  Output 

In  Section  2 of  TR#2  we  developed  formulas  for  the  zonal 
outputs  of  the  third-order  transfer  function.  We  repeat  here  the 
necessary  formulas;  we  are  substituting  the  equivalent  transfer 
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Fig.  2.1  Equivalent  Receiver 


function  K for  wherever  appropriate,  and  making  very  minor 
notational  changes  where  useful. 


P 

X(t)  = total  input  signal  = F._  x (t)  ; 

p X p 

x (t)  = Re [z  (t)  exp  (j2n  v t)  3 ; 
p P P 

z (t)  = z * (t)  ,* 

-P  P 

v = -v  ; 

“?  P 

y (t)  = third-order  output; 

K 


(2-2) 


(2-3) 

(2-4) 

(2-5) 

(2-6) 


yK(t)  - a 


i: 


8 ^3  P ^ ' P 2 * 3 


(t)  exp[j2^(v  + v + v ) t] ; (2-7) 


*Vp2'P3(t) 

W/2 


= rrr  df,df0df,  G (f)  exp[j2nt(f  + f2  +.f3)] 

JJJ  12  3 P1'P2'^>3  1 A . 

-W/2 


z (f.)  z (fJ  z (f,)? 

Pi  1 p2  2 p3  3 


G (f)  = K(f  + v .f,+  v„  .£,+  v ) 

pl-p2'p3 


pl  -2  -p2-'3  p3 


(2-8) 

(2-9) 


The  sum  in  Eq.  (2-7)  contains  8P  terms  altogether,  which 
consist  of  4P3  terms  and  their  conjugates.  Many  of  these  terms 
are  identical  because  they  represent  subscript  permutations. 
Furthermore  most  of  them  will  not  fall  at  carrier  frequencies 
that  can  create  interference;  this  is  the  case  for  any  term  all 

- i .i.  - i ^ A *-  *1  ^rvr  ovnmn  1 o 

01  wnose  bUDbUipua  axe  ^ j.  w y ^ 9 — ‘jr  • — - 
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We  will  now  specialize  Eq. (2-7)  to  those  cases  that  might 
be  of  interest  for  computer  modelling;  we  will  subsequently 
narrow  this  list  down  even  more. 

In  the  tabulation  that  follows  we  will  write  the  outputs 
in  the  form 

y (c)  = Constant  • Re  {a  (t)  exp[  j2TTvmt]  } , (; 

K P1'P2'P3  T 


with 


Px  > 0 

p2  > 0 
p3  < 0 


*3 

H 

i 

'4 


vm  = v + v + v . (2-10< 

T P1  p2  p3 

The  "Real  part  of"  consolidates  terms  in  Eq.  (2-7)  with  their 
conjugates  and  the  constant  takes  into  account  the  number  of 
permutations  that  lead  to  identical  terms  . 

It  will  be  recalled  from  Section  2.4.3  of  TR#2  that  for 
certain  combinations  of  interfering  signals,  many  different 
interactions  will  lead  to  inband  interference.  We  have  decided 
thai  to  keep  the  computer  programs  manageable  it  is  reasonable 
to  require  the  operator  either  to  ascertain  by  inspection  of 
the  relative  power  levels  which  component  is  the  most  significant, 
or  to  run  all  cases  separately  and  combine  the  outputs  afterwards. 

In  all  of  the  cases  of  interest  we  allow  the  possibility 
that  none  of  the  interacting  frequencies  is  the  desired  fre- 
quency. The  computer  program  must  thus  allow  for  the  possibility 
that  the  linearly  amplified  component  .is  distinct  from  any  of 
the  interfering  complex  envelopes. 
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We  now  tabulate  the  cases. 

(2-11) 


(2-12) 


(2-13) 


(2-14) 

2 . 5 Interference  Spectrum 

It  is  useful  to  have  some  idea  of  the  shape  of  the  spectrum 

of  the  interference  envelopes  a (t) . To  this  end,  let  us  write 

£ 

the  time  domain  analog  of  Eq.  (2-8)  ; we  will  drop  some  of  the 
subscripting  where  it  will  cause  no  confusion.  We  have 

a<t)  " JIfdtldt2dt3  9(tl't2'V 

z (t-t  ) z (t-t  ) z (t-t  ) , (2-15a) 

P]  1 p2  2 p3  3 

where  g(tj.,t  , t ) is  any  function  whose  transform  equals  G on 
h Vto  nnKa  r\f  i.on  i.n  E*c[  - * p— 3}  • 


v =v1,v  =v.,v  = - vn 

Pi  1 p2  1 p3  i 


yK(t)  = (3/4)  Ret^  ^ ^(t)  exp ( j2nvT (t)  ] 

v = v ,v  = v ,v  = - v 
P1  1 P2  1 P3  2 

f 

yK(t)  = (3/4)  Refa^  1 _2  (t)  exp(j2nvTt)]  J 

V = V , ,v  = v.,v  = -v„ 

Pi  1 p2  ^ P3  2 

yK(t)  = (3/2)  Re[ax  2 _2  (t)  exp  ( j2TrvTt)  ] 


v - v ,v  = v.,v  = - v. 
Pl  1 Po  2 P-3  3 


yR(t)  = (3/2)  Re[a^  ? _3(t)  exp ( j2TrvTt)  ] J 
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(2-15b) 


„dtldt2at3  9(tl'W 


exp[-j2n(t1£1  + 

= G ( f ^ , f ^ » f when  1^1^  W/2,  i = 1,2,3. 

For  future  reference  we  should  keep  in  mind  that  since  p^<  0, 


z (t)  = zf  I (t)  . 
P3  |P3I 


(2-16) 


Let  R (t)  be  the  autocorrelation  function  of  a(t) 

3. 


R.  (t)  = E[a (t) a* (t+T ) ] . 


(2-17) 


Substitution  of  Eq.  (2-15a)  in  Eq.  (2-17),  with  the  introduction 
of  ne’rf  dummy  variables,  yields 


R 


(f ) J* • . . J*dt^ . . . dt^  ^ (tjL*  t2' t^) 


U(t1> . . . , t^jT) 


42-18a 


where 


U ( t , . . , t ; t ) =E 
1 6 


f z (t-t, ) z*  (t  + t - t ) 


•z  (t-t  ) Z*  (t  + T - t ) 

P2  2 P2  5 

®‘Z  (t-t  ) Z*  (t  + T - t ) , 

L P3  3 p3  6 J 


(2-18b 


As  usual,  by  E()  we  mean  the  "expected  value  of". 


Before  going  any  further,  we  should  observe  that  the  form 

of  Eq. (2-15a)  indicates  that  in  the  case  defined  by  Eqs. (2-11) , 

(2-13)  and  (2-14) , the  interference  is  highly  correlated  with 

tiiG  component  envelope  z \t)  • Consequently,  xn  any  o£  i-liese 

Pi 

three  cases,  if  z (t)  is  the  envelope  of  the  desired  signal, the 

Pi 

power  spectrum  of  “the  interference  may  be  of  little  interest 
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because  the  interference  actually  bears  useful  signal  informa- 
tion. We  will  therefore  exclude  those  cases  from  consideration 
in  evaluating  Eq.  (2~18b) . 

That  equation  is  hopeless  to  evaluate  as  it  stands  because 
Uv-  determination  of  the  expectation  of  the  sixfold  products 
requires  information  on  the  joint  statistics  of  envelopes  at 
six  time  instants.  However,  some  progress  can  be  made  if  we 
assume  that  the  envelopes (zp^  (t) } are  complex  Gauss: an  processes 
with  identical  covariance  functions.  (They  will  be  identical 
processes  when  the  subscripts  coincide.)  We  will  make  this 
assumption,  and  v/rite  the  common  covariance  as 

R (t)  ~ e{z  (t)  z*  (t+T)|,  i = 1,2,3.  (2-19a) 

% C P . P . J 

1 

We  note  that  this  covariance  satisfies 

R (-T)  = R *(t)  . (2-.l9b) 

z z 

It  will  be  convenient  to  approximate  these  autocorrelations 
as  being  associated  with  a rectangular  power  spectrum  of  band- 
width W: 

W/2 

R^t)  = exp  ( j2T,fT)df . (2-20a) 

2 -W/2 

(It  perhaps  should  be  pointed  out  that  we  are  ignoring  the  scale 
factors  on  these  autocorrelation  functions  and  power  spectra, 
which  are  immaterial  to  the  shape  of  the  interference  spectrum.) 

It  will  be  observed,  that  where  convenient  we  may  use  equally  well 

I*  W/2 

R (t)  = i exp  (- j 2p£t ) d£  (2-20b'/ 

2 -W/2 

because  the  R (t)  defined  by  Eq.  (2-20a)  is  pure  real, 
z 
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In  the  appendix  to  this  section  we  derive  the  general 
form  of  expectations  of  the  type  in  Eq.  (2-l8b) . The  results 
these  may  be  used  in  conjunction  with  Eq.  (2-18)  to  derive  the 
autocorrelation  function  of  the  interference  envelope,  and  hence 
the  spectrum.  We  wish  to  restrict  attention  here  to  the 

special  cases  described  by  Eqs. (2-12)  and  (2-14). 

In  Eq. (2-12)  there  are  two  distinct  carrier  frequencies 
and  v^,  and  consequently  two  distinct  interferer  envelopes 
so  that 


z (t)  = z (t)  = z , (t) 


1- 


z (t)  = z*  (t) 
p3  2 


Equation  (2-18b)  then  becomes 


U(t_  t • • * / 1 ..7  T ) 
1 6 


J 


= E 


7i(t-ti)  vfc-y  *2(t+T-t6) 


^Z1*(t+T~t4)  Z1*(t+r-t5)  Z2*(t-t  )j 

which,  from  Eq.  12-A7)  in  the  appendix  is 


U(t^,..», tg } T ) 


R (t,-t  -t) 
Z 6 3 


vt-v  v vt-v  v 


(2-21) 


(2-22) 


(2-23 


+R2(T-t4+  t2)  R2(t-V  t]L) 

If  we  now  assume  the  autocorrelation  function  given  by  Eq. (2-20) , 
this  may  be  substituted  in  Eq,  (2~18a)  to  yield 
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W/2 

Ra(7)=  III  dVf2df3  I Jat! — at&  S'WV  9*(t4't5't6) 

- W/2 

j'eXpr-j2-t(f1t1-£1t4+f1T+i2t2-f2t5+f2-r£3t3-f3t6+f3T)], 

|+exp[-j2r(£1t1-£:!t5+f1T+(2t2-£2t4+£2T+f3t3-f3t6+f3T)]| 


Evaluation  of  the  integrals  in  t^,...,t^  yields 
W/2 

Ra(T)=  JJJ  df^df^df^  exp[-i2TTT  (f^  + f 2 + f^'] 
-W/2 

G(fl,£2.£3)  G*  (£rt2.f3)' 


+G(£1.£2,£2)  G*  (f2.f1.f3)J 


Since 


G(f2'£1’f3)  = K(,2+Vfi+vrVv2) 


Mf1+v1.£2+v1.f3-v2) 


(2-24) 


(2-25) 


= G(frf2,f3)  , 

the  large  bracket  in  Eq„(2-25)  is  actually  equal  to  twice  the 
first  summand  in  it.  We  can  v/rite  down  by  inspection  the  power 
spectrum  of  the  interference  as 

CD 

G (f)  = I dr  exp(-j2TTfT)  R (t ) 
a I®  a 

W/'2 

- 2 LIT  dfl  df2  af3  6'*-VVf3> 

-W/2 

I fZ  f -F  f \ I ^ 

1 — i'  ~2'  ~y  1 
or 


2-11 


where 


Vf)“  2 JJ  dfiaf2l<Mfl'f2<f-frf2>*2 

(fyf2)pMf 

nf  = f(fX'f2):*fiJ  * w/2'!f~Vf2'  * W2)'  (2‘ 


w-hen 


v » V = V V = -V 

P1  P2  1 P3 


(2- 


We  are  particularly  interested  in  the  tail  of  the  spectrum 


wnere 


W/2  * I f | * Tv'2, 

in  comparison  with  the  peak  value.  Suppose  we  let 


G - Min  | G (f . , f _ , f _ ) 

m | f . i £W/2  1 2 3 


(2- 


and 


G = Max  i G (f , , f f ) 
M |f.|*W/2  1 2 3 


(2- 


Then,  for  any  f 


G (f)  * 2 G 
ei  M 


dfldf2' 


(2- 


G (f)  * 2 G 2 f!  df.  df_. 
a m * «L  1 2 


a. 


(2- 


The  integral  over  0^  can  be  evaluated  fairly  readily;  we  find 


*■  ?'  2 
( 3W  /4  ~ f 


* If | £ W/2 


p i ' 

I df.df  * \ . , 2 , > 

ft  “ 2 (J3W/2  - j f | ) /2  , W/2  <|f|^  3W/2J 


(2- 


-26b) 

-26c) 

-27a 

-27b 

-28a 

-28b 

29) 
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From  Eqs(2-28)  and  (2-29)  we  can  bound  the  ratio  of  inter- 
ference power  in  the  channels  centered  on  v + W to  the  inter- 
ference power  in  the  channel  centered  on  v . We  have 


adjacent  channel  interference  power  M 
direct  channel  interference  power  4G 

m 


(2-30) 


We  can  repeat  the  whole  procedure  now  for  the  case  in  Eq. 
(2-14)  where  all  three  carriers  are  distinct: 


v = v ,v  = v , v. 

Pi  1 Po  2 I 


'3* 


(2-31) 


The  expectation  U is  given  by 


U(t  , . . .,t  ?t) 
1 6 


= E 


'*l(t'V  *2(t'V  Z3(t+T'V 

z,  (t+T-t.)  Z_  ( t+T  — t ) Z _ ( t— t ) 


(2-32) 


from  which 

U(t. t ;T ) = R(T-fc+t.)  R(T-t_  + t_)  R(t-t-r).  (2-33) 

J b 41  5263 

A comparison  with  Eq.  (2-23)  shows  that  the  interference  spectrum 

is  exactly  one-half  that  found  in  the  previous  case  so  that  Eq. 

(2-30)  applies  in  this  case  also. 


2.6  Final  Selection  of  Interference  Combinations 


We  are  now  in  a position  to  make  a selection  of  the  inter- 
ference cases  to  be  modelled.  We  recall  that  Eqs.  (2-11)  to 
(2-14)  defined  four  basic  combinations  of  interfering  frequencies, 
and  for  each  of  these  cases  it  would  in  general  be  possible  to 

examine  situations  where  the  interference  carrier  v;  fell  on 

T 

either  the  desired  carrier  v or  the  adjacent  channel  carriers  v+W. 
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The  case  given  by  Eq.  (2-11)  involves  interaction  of  a 
signal  with  itself  so  that  the  interference  band  is  centered 
on  the  same  carrier  frequency.  The  only  two  situations  of 
interest  are  those  where  = v and  v^=  v + W.  (The  lower  ad- 
jacent channel  case  i^  essentially  identical  to  the  upper  ad- 
jacent channel  and  need  not  be  treated  separately.)  We  find 
it  reasonable  to  assume  that  a well-designed  AGC  circuit  will 
preclude  significant  self-interference  of  the  desired  signal. 
We  therefore  restrict  this  case  to 


Case  1: 


vm  = v + W . 
T 


(2-34) 


The  situation  in  Eq.  (2-12;  is  that  of  carriers  at  and 
producing  an  intermodulation  carrier  at  2v^  ~ v 2°  We  a^a^n 
assume  that  an  adequate  AGC  makes  the  case  = v uninteresting. 
This  leaves  only  the  question  of  whether  to  allow  modelling  of 
the  cases 


- 2v^  - v + V7, 


as  well  as 

v?  = 2v . -v. 

The  conclusion  we  draw  from  Eq.  (2-30)  is  that  the  adjacent 
channel  interference  effect  is  unlikely  to  be  as  strong  as  the 
direct  channel  interference  effect  so  that  for  any  reasonably 
well  behaved  K we  can  restrict  attention  to  the  situation  where 
the  intermodulation  carrier  falls  on  the  desired  carrier.  We 
are  thus  led  to  take  as  the  second  case: 


Case  2:  v = v = v, ;v  = v-2v. ;v  = v.  (2-J5) 

P1  p2  1 P3  1 T 

The  situation  in  Eq.  (2-13)  is  a cross-modulation  inter- 
ference where  the  interference  carrier  is  at  the  frequency  of 
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one  of  the  two  interfering  carriers.  The  classic  case  in- 
volves crossmodulation  of  the  desired  signal,  which  certainly 
must  be  evaluated.  However,  it  woulo  appear  equally  important 
to  consider  the  situation  where  a moderately  strong  adjacent 
channel  signal  is  splattered  into  the  desired  band  as  a result 
of  crossmodulation  by  a strong  out-of-band  signal.  We  are 
thus  led  to  two  more  cases: 


Case  3:  v 

pl 


(2-36) 


Case  4:  v 

pl 


V_  = V+W,  V 
T P2 


(2-37) 


The  final  situation  is  that  of  Eq. (2-14)  involving  three- 
frequency  intermodulation.  We  note  that  Case  2 and  Case  4 above 
are  special  cases  of  Eq. (2-14) . Furthermore,  from  the  results 
of  Section  6 of  TR#3  we  know  that  the  probability  that  three 
carriers  have  sufficient  power  to  produce  this  interference 
is  small  compared  to  the  probability  of  the  two-carrier  inter- 
actions described  by  Eqs. (2-35)  and  (2-37) . We  therefore  ex- 
clude this  option  from  the  modelling  capability. 

In  summary,  the  cases  described  by  Eqs.  (2-34)  through  (2-37) 
will  form  the  basis  of  our  modelling. 

2.7  Scaling 

The  complex  envelopes  appearing  in  the  several  expressions 
include  implicitly  a scale  factor  proportional  to  the  square 
root  of  the  nominal  carrier  power.  These  implicit  scale  factors 
are  multiplied  together  (also  implicitly)  in  determining  the 
peak  voltage  of  the  complex  envelope  of  the  interference. 
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Considerable  efficiency  can  be  achieved  in  the  operation 
of  the  computer  programs  if  these  scale  factors  are  made  ex- 
plicit and  their  product  used  to  scale  the  interference  output 
after  it  has  been  calculated.  In  this  way  the  effect  of  varying 
power  level  can  be  determined  by  scaling  a single  output  sequence 
rather  than  by  repeating  the  entire  triple  summation  involved 
in  the  tapped-delay  line  model. 

It  is  also  appropriate  to  include  the  peak  magnitude  of 
K^l'^2'^3^  fc^s  scaling  so  that  the  tapped-line  co- 

efficients have  a relatively  restricted  set  of  magnitudes. 

2 . ft  Tap  Spacing  and  Sample  Spacing 

We  would  like  to  discuss  here  the  interrelation  between 
the  tap  spacing  and  the  sample  spacing.  Let 

6 = time  interval  between  adjacent  samples  of  the 

complex  envelopes;  (2-38) 

t - time  interval  between  taps  in  tapped  delay  line 
° model.  (2-39) 

The  range  of  values  of  t is  determined  by  the  requirement  of 

getting  a good  fit  to  the  transfer  function.  Values  of  t near 

zero  will  be  used,  for  example,  in  approximating  polynomial  fits 

to  the  transfer  function.  At  the  other  extreme  we  can  assume 

that  t < 1/W  which  is  the  largest  value  that  permits  a Fourier 
o 

representation  of  Kff^f^f^)*  We  thus  have 

0 < t < 1/W.  (2-40) 

o 

A second  constraint  is  imposed  by  compatibility  with  the 
sampled-data  representation  of  the  complex  envelopes.  We  clearly 
must  constrain  t to  be  an  integer  multiple  of  6: 

t = integer* 6.  (2-41) 

o 
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The  sampling  interval  itself  is  constrained  by  the  need 
for  adequate  representation  of  the  complex  envelope  of  the 
interference.  Since,  by  the  earlier  assumptions,  this  envelope 
has  a spectrum  occupying  the  interval  (-3W/2,  3W/2)  a sampling 
rate  of  3W  samples/second  is  the  minimum  allowable  to  permit 
Nyquist  sampling.  To  avoid  the  need  for  (sin  x./x)  sampled 
data  filters  it  is  wise  to  allow  at  least  some  margin  and  re- 
quire that 

6 < 1/4W.  (2-42) 


For  this  initial  effort,  which  must  be  looked  on  as  a 
validation  of  the  possibility  of  computer  simulation  of  the 
models,  we  have  chosen  to  satisfy  all  of  the  constraints 
simultaneously  by  requiring  that 


6 = t 

o 

and 


(2-43) 


0 < t < 1/4  W.  (2-44) 

o 

The  penalty  imposed  by  this  lack  of  flexibility  in  setting  6=t 
is  a restriction  of  the  "wildness"  of  the  transfer  functions 
that  can  be  accommodated  by  the  model.  We  should  point  out, 
however,  that  this  same  restriction  greatly  reduces  the  measure- 
ment or  computational  burden  in  determining  the  values  of  the 
transfer  function  on  a cubic  lattice. 
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APPENDIX  TO  SECTION  2: 

CALCULATION  OF  TRIPLE  MOMENT  OF  COMPLEX 
GAUSSIAN  RANDOM  VARIABLES 


Let  x^, . 
the  moment: 


xe  be  complex  Gaussian  variables  and  let  M be 


R = E(Xlx2x3x*x*x*). 


(2-A1) 


By  a Gram-Schmidt  procedure  we  can  represent  {x  J as  a trans- 
formation  on  uncorrelated  unit  variance  variables  in  the  form 


where 


d13. 

a2ly'+a22y2 


6 * a61yl  + 


+a66y6 


6 , the  Kronecker  delta, 

mn 


(2-A2) 


(2-A3) 


If  we  substitute  Eq.  (2-A2)  in  Eq.  (2-Al) , the  result  is  a sum 
containing  6!  products  of  the  jy^j  and  their  conjugates.  It 
would  be  bad  judgment  to  write  this  sum  out  because  most  of  the 
terms  vanish  when  we  take  the  expectation.  In  fact  the  only  non- 
zero expectations  are  those  of  the  form 


in  which  i,j,k  may  or  may  not  be  distinct  subscripts.  We  see 
immediately  that  Y^'V^  and  Y&  will  never  enter  into  the  calcula- 
tion of  R because  they  never  appear  unconjugated  in  the  product. 
We  will  need  tne  three  moments  of  the  v nit-mean  exponential 
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distribution  in  writing  down  tbe  expected  products: 

E{ 1 yi l 2p}  = pJ- 

We  then  have 

R = 6(aua21a3i)(a41a5ia61)* 

+2 (alla22a31+alla21a32) (a51a61a42+a41a61a52+a41a51aS2» ‘ 

+2  (aUa2la33)  (a51a61a43+a41a61a53+a41a51a63)  * 

+2  (aUa22a32)  (a51a62a41+a42a62a51+a42a52a61) 

+ (alla22a33)  <a4laS2a63+a51a62a43+a61a42a53+a41a53a62 

+a51a63a42+a61a43a32) 

This  latter  expression  needs  to  be  rephrased  in  terms  of 
the  covariances  of  pairs  of  jx  j: 


R = E-ix  x > m 7.  a a . 
mn  l m n J i mi  ni 


Now  if  X^,  X^^and  X^  were  independent  and  if  X^,  Xr , X^  were 
permutation  of  them,  then  R would  contain  a product  of  the  form 

R,  „ R0  R-, 

1,P1  2 ' P2  3'p3 

where  (P^P^P^  was  a permutation  of  (4,5,6).  We  are  thus  led 
to  conjecture  that 

R = R14  R25  R36  + R15  R26  R34  + R15  R24  R35 

+ R14  R26  R35  + R15  R24  R36  + R16  R25  R34*  <2-A7) 

If  we  write  out  this  sum  using  Eq.  (2-A6)  we  can  verify  its 
coincidence  with  Eq.  (2-A5)  and  thus  correctness  of  our  conjecture. 
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SECTION  3 


DETERMINATION  OF  MODEL  PARAMETERS 


We  start  here  with  Eq.  (2-10)  of  the  last  section  which  we 
repeat  here: 


yK(t)  = constant  -Re{a  > (t)  exprj2rrv  t]  ] 


P1'P2'P3 


where 


P1  > °'\ 

p2  > ° ; 
p3  < o 


V = V + V + V 
T P1  p2  ?3 


The  complex  envelope  a (t)  in  Eq.  (3-la)  is  given  by  Eq.  (2-8) 

jeL 

as 

W/2 

a£(t)  = JJJ  d-  G£(^}  exprj2nt(f1+  f2  + f^)] 

-W/2 


•ZP1(fl)  Zp2(f2)  Zp3(f3}' 


in  which 


then 


G (I)  = K(f  + v , f + v , f + v ). 
fi  1 Px  2 p2  3 p3 

Our  modelling  procedure  relies  on  the  property  that,  if 

CM!)  « G(f)  for  |f.j  < w/2,  i * 1,2,3, 

CD 

a£(t)  * JJ7d-  'S(I)expCj2^t(f1+  f2+  f2)] 

m 

\{fl]  Zp2(f2>  Zp3(£3>- 


(3-4) 


3-1 


by  virtue  of  the  bounded  support  of  the  input  spectra  [z  ).  We 

^i  . 

have  furthermore  chosen  to  restrict  attention  to  approximations 
6(f)  of  the  form 

6(f)  = „v  nv  J B(n  ,n  ,n  )expT~j2n6(n  f + n„f  + n f )]. 

i z j 11  £.  £ jo 


Substitution  of  Eq~  (3-6)  in  Eq.  (3-5)  then  leads  to 


(3-6) 


a ( t ) = E y E B(n  ,n  ,n  ) 
E n!  n2  n3  123 


3 ■ 00 

. T1,  f df.  exp  rj2rrf.  (t-n.6)]  z (f.) 
II  J i i i P-  i 

i=l  -°° 


(3-7) 


or 


a (t)  = E T E B (n  ,n  »n  ) Tj  z (t-n.6). 
E n!  n2  n3  1 2 3 p.  i 


(3-8) 


3.1  Choice  of  Model  Parameters:  General  Considerations 


Referring  to  Eq. (3-6)  the  parameters  that  must  be  specified 


are 

a)  the  tap  spacing  6 

b)  the  range  of  indices  {n^n^n^}  in  the  sum 

c)  the  coefficient  set  [B  (n^n^n^)  } . 

Fox  computational  purposes  it  is  reasonable  to  require  that  the 
set  of  indices  be  identical  for  each  coordinate  so  that  Eq. (3-6) 
can  be  specialized  to  read 


6(f) 


TIT. 

rn-p 


N + 1 5 


n. 


SNI+ 


LV 


B(Vn2 

1 


e 


- j2Tr&n . f . 

li 


i =1 


(3-9) 


3-2 


The  analogous  time  response  of  Eq.  (3-8)  then  becomes 


a (t)  = HF 

e «-  nr 


N +1  * 


n. 


U"3J 


B<ni-n2'n3)  n v(t-ni6) 

i=l  1 


(3-10) 


£ N + N 


In  these  expressions  N is  then  the  number  of  delays  used  in  the 

model  for  each  input  sequence,  while  the  number  can  be  thought 

of  as  an  overall  delay  of  the  output  sample  sequence  relative  to 

the  linear  components  of  the  model.  The  number  N can  be  assumed 

to  be  fixed  ahead  of  time  by  complexity  limitations  on  the  com- 

3 

puter  programming.  It  is  therefore  necessary  to  determine  N + 2 
model  parameters: 

ci)  the  delay  spacing  6 

b)  the  overall  delay  N 

3 ■** 

c)  the  N values  of  the  {b  (n^,n0,n^) } . 

If  the  function  G(fj  were  specified  at  all  values  of  f by  an 
analytic  description,  it  would  be  reasonbly  straight forward  con- 
ceptually to  find  a choice  of  the  N^+  2 parameters  that  minimized 
the  quadratic  approximation  error 


JJJ  df|G(f)  - G ( f ) | 2 . 

What  we  will  actually  be  working  with,  however,  is  a finite  set 
of  calculated  values  as  measurements  of  G(_f)  in  the  form 


G(m.  f , m f , m f ) 

1 o 2 o 3 o 

as  m^.m^,  and  m^  range  over  some  small  set  of  integers.  If 
the  number  of  measurements  per  frequency  coordinate  so  that 

3 

total  number  of  measurements  is  M , then  clearly  we  require 


M is 

fho 

that 


3-3 


exceed  the  number  of 


so  that  the  number  of  measurements  wi LI 
parameters  to  be  determined.  Thus  we  must  take 


M * N + 1. 


3 

On  the  other  hand,  the  determination  of  each  of  the  M data  points 
will  require  either  a significant  measurement  effect  or  significant 
computation  time.  Therefore  we  believe  that  it  is  reasonable  to 
use  only  the  minimum  number  and  hence  require  that 


M = N + 1, 


(3-lla) 


where  N is  the  largest  value  of  N for  which  modelling  is  to  be  done. 


At  this  point  we  should  observe  that  the  smallest  "interesting 
value  of  N is  N = 2 corresponding  to  a two- tap  model,  or,  in  the 
limit  of  small  tap-spacing  6,  a frequency  power  series  with  linear 
terms  in  each  frequency.  We  have  chosen  to  allow  values  of  N as 
large  as  4 which  provides  considerable  flexibility  beyond  the 
minimally  interesting  model;  the  corresponding  value  of  M = 5 

which  requires  125  measurements  or  calculations  is  probably  as 
large  as  can  be  conveniently  accomplished  with  any  reasonable 
economy  of  either  computer  or  measurement  time.  We  thus  will  re- 
strict attention  to 

M * 5 (3-12a) 


and,  from  Eq.  (3-11) , 

N s 4.  (3-12b) 

The  most  favorable  location  of  the  frequency  lattice  points 
is  not  immediately  apparent.  To  be  more  specific,  it  is  not  ap- 
parent how  close  to  the  band  edges  of  + W/2  the  extreme  data  points 
should  be.  Although  an  argument  could  be  made  that  choosinq  the 
frequencies  of  + W/2  as  2 of  the  data  coordinates  on  each  dimen- 
sion gives  undue  weight  to  possibly  anomalous  band-edge  phenomena. 


'v  K,swe^' ^fr**1*^ 


this  choice  is  a conservative  one  in  that  it  will  make  any  model 
weaknesses  most  apparent.  We  therefore  will  henceforth  assume  that 
the  measured  or  calculated  data  consist  of  the  values 


for 


Y (m^m  #m  ) 


r-  w/2  + 


G 


-W/2  + 


L -W/2  + 


(m^DW/N,-. 
(m2-l)W/N,  , 

(m  -1)W/N  . 


1 £ m^,m2,m2  £ N+l. 


(3-13a) 


(3-13b) 


We  now  define  Y in  the  obvious  way  as  the  value  of  the  ap- 
proximating function  G at  the  same  lattice  frequencies,  and  define 
an  error  criterion 


N+l  N+l  N+l 

V=  E E T.  |Y(m1,m2,m3)  - YOn^m  ,m  )|  (3-14) 

m =1  m =1  m =1 
12  3 

v;hich  is  a discrete  version  of  the  quadratic  error  criterion.  It 
should  be  observed  that  V is  implicitly  a function  of  the  param- 
eters N , 6,  and  {b (n^/n^.n^) ] . The  objective  of  the  parameter- 
extraction  program  is  then  to  minimize  V by  the  choice  of  these 
implicit  arguments. 

Before  proceeding  further  it  is  helpful  notationally  to  let 
P(n  ,n  ,n  ) = B(n  +NJ#n  +NJ#n  +1^)  (3-15) 

so  that  Eqs.  (3-9)  and  (3-10)  can  be  rewritten  as 
G(f)=  exp  C-j2TTNI&  (f1+f2+f;J)  ] 


N 


N 


N 


^ a / 


H VII 


V1  n2=1  n3=1 


1/U2 


'V 


3 

0T“t 

I i 

i=l 


- j2n6n . f . 

J ii 


l3-ibj 


3-5 


(3-17) 


a (t)  « 
£ 


N 

E 


^ ^ P (n  ,n  ,n  ) f*  2 (t~n. 6-N_6) . 

»,-l  n =1  1 2 3 Pi  1 * 


It  can  be  verified  that  for  any  fixed  6 and  N„,  the  determina- 
tion  of  the  best  values  of  ^(n^,n2,n3^  •*'s  a routine  quadratic  mini- 
mization. (We  will  present  the  derivation  of  these  coefficients 
subsequently.)  The  problem  thus  reduces  to  finding  an  efficient 
algorithm  for  determining  the  best  choice  of  6 and  N 

X 

3 . 2 Dependence  of  v on  6 

Although  it  is  theoretically  possible  to  find  the  optimum 
value  of  6 by  differentiating  the  quadratic  error  V with  respect 
to  6 and  equating  this  derivative  to  zero,  it  appears  to  be  more 
realistic  computationally  to  search  for  the  minimum  V by  evaluat- 
ing Eq. (3-14)  for  several  discrete  values  of  6.  This  approach 
also  guarantees  that  the  value  of  6 we  select  will  be  an  approxi- 
mation to  the  value  yielding  a global  minimum  rather  than  one 
which  yields  only  a local  minimum.  We  now  need  to  investigate 
the  limits  of  this  computer  search,  or,  more  exactly,  the  lower 
limit,  since  we  have  already  determined  that  6 < l/4w.  We  thus 
need  to  investigate  the  behavior  of  the  optimizing  solutions  as 
6-0  and  then  determine  the  largest  possible  positive  6 which  per- 
mits approximating  this  limiting  behavior. 

From  Eqs.  (3-13)  and  (3-26)  ve  have 

V (m1,m2,m;i)=  exp{-j2-Nl6w[-3/2-3/N+(m1+m2+m3)/N]  } 

N N N 

* E £ £ 0(n  ,n  ,n  ) 

"r1  v1  vl 

rj 

yj> 

exp[-j2TT&niW(-l/2  -1/Nvnu/N)  ] . (3-18) 

i=l 


3-6 


We  clearly  must,  have  each  of  the  Y approaching  a limiting  value 
as  6 - 0;  we  will  denote  this  limit  by  the  subscript  zero: 

Y_(m  ,m  ,m  ) = lim  Y(m.,,m^,m  ).  (3~19) 

0123  , _ 1 2 3 

o-O 


In  Eq.  (3-18)  , let  us  denote  the  triple  sum  by  a (m^.m^/m^)  .. 
By  writing  the  product  of  exponentials  as  the  exponential  of  the 
sum  of  arguments,  o can  be  written  as 

N N N 

o (m  ,m  , m ) = E v £ p(n  ,n  ,n_)  • 

1 Z J m - - Iz3 

V1  n2=1  n3=1 

3 

.exp!~-j2TT5w  E n.  (-l/2-l/N+m./N)1.  (3-20) 
*-  i=l  1 ~ 1 “ J 


If  we  now  expand  the  exponential  in  a power  series  in  6,  the  sum- 
mation in  the  power  series  can  be  commuted  with  the  triple  finite 
sum  on  {ni'n2'n3^  to 

00 

o(m  ,m  ,m  ) = E (-j2nw6) 
k=0 

N N N 

* E T.  E P(n1*n2'n3)  * 

"l*1  n2=1  n3=1 


\\ 

1_1=1 


n . 

i 


m. 


. _,k 

N / J * 


(3-21) 


Now,  if  the  coefficients  O (n^,n^,n^) } were  constants  independent 
of  6,  then,  for  sufficiently  small  6,  the  leading  term  in  the  in- 
finite sum  would  be  the  dominant  term,  and  a (m^m^m^)  would  be 
independent  of  m^m^,  and  m^.  However,  we  can  expect  that  the 
dependence  of  the  {P (nn ,n^,n^) } on  6 will  be  reflected  in  rela- 

4.  4.  ^ 

tions  of  the  form 


3-7 


£ (3  (n  ,n  ,n  ) -*  constant  as  6-0, 

WV 

T.  n.  ft  (n  , n ,n  ) - CQnstant  as  6-0( 

. x 1 1 2 3 0 

ni' n2'  r^) 

r>  o / v constant  t ^ 

£ n.n.  p (n  ,n  ,r.  ) - as  6-0, 

/ \ ^ D 12  3 1 2 

(n1,n2,n3)  6 

and  so  forth.  It  will  thus  be  possible  to  obtain  polynomials 
in  m.,m0,and  m as  6-0. 

J.  Z,  J 

Suppose,  then,  that  we  postulate  that  each  3 has  a finite 
Laurent  development  of  the  form 

P(ii  ,n  ,n  ) = pQ  (r?1-n  , n^)  + p (n1,n2,n3)  6_1 


+. . .+  Pj  (n  ,n2#n  ) 6 


-L 


(3-22) 


+ terms  of  the  order  of  6 or  smaller. 


If  we  let 


PD^ni'n2'n3^  = homogeneous  polynomial 

^ J in  (n^,n2#n3)  of  degree  D,  (3-23) 


then  clearly  we  must  have 


z MVn  .n  ) PD(n  .n  ,n  ) - 0 if  D < l 
(n1,n2,n3)  (3-24) 


for  otherwise  there  would  be  infinite  values  of  a as  6-*0.  Equa- 
tion (3-24)  imposes 


^•6+2^  _ t (~6~f  1)  (-6+2) 


linear  constraints  on  the  [ (n^ , n2 , n3 ) ) corresponding  to  the 
number  of  types  of  homogeneous  polynomials  of  degree  < t.  Since 

3 

there  are  only  N cf  these  coefficients,  it  follows  that 


3-8 


(3-25) 


L (L+l)  (L+2)  £ 6N3 

L S 1 if  N = 1 

L ^ 2 if  N = 2 

L ^ 4 if  N = 3 

L ^ 6 if  K = 4 J 


(3-26) 


The  limiting  value  of  a(  ),  which  we  will  denote  by  o ( ) , is  then 

L i ° 

o (m  ,m  ,m  ) = E (-j2"W) 

O -L  ^ J j_  q 


N 

E 


N 

E 


N 


E Pi(n1,n2,n3) 
nl=1  V1  n3=1 

3 . ..  m. . _-t 

["  ni(-  i - i + f )]  • <3-27) 

i=] 


The  terms  in  this  expression  can  be  rearranged  to  yield 

al  n2  °3 

o (nt,  .m  .m„)  = Y.  c m,  m„  m„ 

o 1 2 3 aeA  a 1 2 3 


(3-28a) 


where 

a = (a1*a2'a3^  (3-28b) 

and 

A = (a:  Ea.  ^ l}.  (3-28c) 

i 

The  number  of  coefficients  in  Eq.  (3-28a) , or,  equivalently, 

the  cardinality  of  A,  is 

/L+3\ 

Card  (A)  = ( 3 J.  (3-29) 

This  might  suggest  the  inference  that  there  are  more  degrees  of 
freedom  in  the  limiting  case  than  in  any  of  the  expressions  for 
nonzero  6.  This  inference  would  be  false,  however.  In  addition 
to  the  linear  constraints  of  Eq. (3-24)  there  are  linear  constraints 


3-9 


™ * * ^ **  **  ■*  ft* « f - j:  ^ jjrf7?j*^»  'i  rivF^.W,''-, , ^ 


vv*-*?y»*-> ^*v v^ 


between  the  different  {P^}  and  the  measurement  valies.  In  fact, 
if  6^  is  some  value  of  6 sufficiently  small  so  that 


a,  (m  , rn  , m ) 
61  1 2 3 


o (m  ,m  ,m  ) 
0 12  3 


3 

then  the  fact  that  there  are  only  N independent  coefficients  in 

3 3 

Eq.  (3-20)  implies  that  of  the  M = (N  +1)  values  of  05  , only 

3 1 

N are  linearly  independent. 


We  now  need  to  investigate  the  race  of  approach  to  the  limit- 
ing form  of  Eq.  (3-28a) . To  this  end,  let  us  define 

L -t  r 

r6 {1)  = n -tio  h {-]  6 exp[-j2TT6Zn.f.J.  (3-30) 

The  rate  of  approach  of  r^(f)  to  its  limiting  value  is  at  least 
as  good  as  would  be  obtained  if  the  summation  on  t included  positive 
powers  in  6 . Wa  find 


lim  6 ^T.  p (n)expf-j2n6En.  f .] 

C rt  t xi 

6-0  n 


ptd) 


(3-31) 


where  is  some  homogeneous  polynomial  of  degree  t in  3 variables. 
It  is  actually  one  part  of  the  frequency  polynomial  approximation: 


v 

L 

l 


vV(-j2n) 

12  3 


tl  l 

„ 0 - . 1 2 3 

E Pt<S)nl  n2  n3  ; 


(3-32a) 


in  which 


(3~32b) 


It  is  not  difficult  to  verify  that  the  slowest  convergence  is 
obtained  when  l = L.  We  will  therefore  restrict  attention  to 
that  case. 


3-10 


w 


rV  r3iWKI^^Wd>>*'5^V*«  i ',tH*,‘  '*-  '**"  * " * ‘;'<,v  ’ ' ' 


We  now  observe  that 


m,  m„  in 
12  3 

n,  n„  n„  = 
12  3 


(2^j) 


bill 


‘Vl^Wa1”  , „ „ 

1+m^  l+m^  1+irij  Z1  Z2  Z3 

Z1  2 2 Z3 


(3-33) 

where  all  three  integrals  are  on  the  unit  circle.  Hence  any 
homogenous  polynomial  of  degree  D in  n can  be  approximated  by  a 
linear  combination  of  polynomials  of  the  form  (n^z^+n^^n^  )D. 
The  linear  constraints  of  Eq.  (3-24) can  therefore  be  phrased  as 


I 


1 


2 P (n)  (n  z +n  z +n  z ) = 0 
n L 112233 

for  every  nonzero  2:  and  for  every 
D < L. 

Suppose  chat  we  define 

9(w,2)  = E 57  >:PL(n)(ni*1+n2z2+n3Z3)D- 
D=0  ii 

Transposing  the  order  of  the  two  sums, 

g(w,z)  = ^PL(n)exp  Iw^z  +n2z2+n  z3)]  . 


(3-34) 


(3-35) 


(3-3<\ 


The  constraint  equation  now  becomes,  from  Eqs. (3-34)  and  (3-35) 
-L 


Lim  w g(w,zj  < ® for  every  nonzero  z. 
w-0 


(3-37) 


We  observe  now  that  the  right  hand  side  of  Eq. (3-36)  is  a 

polynomial  in 

wz^  wz2  wz3 

e , e , and  e 

This  polynomial  is  of  degree  N in  each  variable  with  t-h<?  zero 
degree  terms  in  each  variable  missing.  It  is  readily  seen  that 
the  only  polynomial  of  this  form  that  satisfies  Eq.  (3-37)  is 


3-11 


L 

g(w,z)=  exp(w  S Z .)  ft  <3-38) 

m=l 

up  to  a constant  multiplier.  In  this  expression  the  [p  are  non- 
negative integers  with  at  least  one  p ^ nonzero  for  each  m.  Further 
more  since  the  product  includes  a term  of  the  form 

* . — i 

exp  ! w E E P . z . ! , 

L.  mi  mi  ij 

the  exponent  in  this  product  must  he  no  greater  than 
(N— 1 ) (z^+  z2+  z^) . 

W'3  therefore  have 


E p . >0  for  every  m. 

mi 

i 


(3-39a) 


E P . ^ N-l  for  every  i.  (3-39b) 

m mi 

The  coefficients  (9  (n) } could  be  found  by  writing  out  the 

Li 

L-fold  produce  in  Eq. (3-38)  and  matching  coefficients  with  Eq. ( 3— 3f 
However,  it  is  much  more  to  the  point  to  compare  Eq. (3-30)  with 
Eq.(3-36)  for  the  special  case  where  0^ (n)  is  zero  for  l < L,  and 
note  that  in  this  case 

r6(f)  = &"Lg(-j2TT6,f)  (3-40) 

so  that  from  Eq. (3-38) 

r&(f)  = & L exp(-j2n6Ffi) 

L 

Fl  {l-exp[-j2TT6rp  .f.]}.  (3-41) 

m=l  i 

The  limiting  form  of  this  can  be  written  by  inspection  as 

L 3 

r (f)  = (j2tr)L  r E p .f..  (3-42) 

o I L . - mi  l 

m=l  i=l 
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Let  us  write 


cp(6)  = r6(f)/ro(f) 

in  order  to  evaluate  the  approach  to  the  limit.  Then 


(3-43) 


cp  ( 6 ) = exp  (-j2Tt6v'f  . ) 


L 1-exp [-j2TT6r  u .f.] 

ft  - 1 mi  • 


(3-44) 


m=l 


j2T'&IH  . f . 
l mi  l 


The  multiplier 


exp  (- j2TT6Ef^) 

can  actually  be  subsumed  in  the  multiplier 
exp(-j2TTN,6  £f  } 

which  appears  in  its  discrete  version  as  the  first  multiplier 
in  (Eq.3-18).  We  can  therefore  take  the  approach  to  the  limit 
to  be  that  of 


_ 1-exp [-j2n6L’p  .f.] 
L ^ J . mi  i 


<p(6)  = Jj 
m=l 


j2n6”u  .f. 

. mi  i 


(3-45) 


which,  for  sufficiently  small  6 is  approximately 

$(6)  « 1 - jnft  r l U . f . . 

J m i mi  l 


(3-46) 


The  absolute  fractional  difference  between  Fc (f)  and  T (f)  can 

o — o — 

then  be  taken  to  be 


e(&,f)  = n&  T,  r h . f.  . 

— in  i mi  l 

This  is  clearly  largest  at  the  band  edges  so  that 


el  £ (TT6w/2)  ~ 


u 


m i 


mi 


(3-47) 


/ *5  A Ci\ 


From  Ea. (3-39b)  this  can  in  turn  be  bounded  above  by 


1 


3-13 


|e|  s 3{K-1)tt6w/2, 


(3-49) 


or,  using 


H s N = 4, 


(3-50) 


{e|  * 9n6w/2 , 


(3-51) 


It  is  a numerical  convenience  to  have  the  smallest  value  of 
6w  be  an  integer  power  of  1/2.  At  6w  = 1/128  we  have 

U|  < 9tt/256  if  6w  = 1/128.  (3-52) 

This  represents  about  a 6°  discrepancy  from  the  limiting  value 
when  this  discrepancy  is  90°  out  of  phase  as  suggested  by  the 
form  of  Sq.  (3-46) . 


We  can  also  check  the  relative  amplitudes  using  Eq. (3-44) . 


We  have 


_ | sin  n&  T.  p . f . | 
L mi  l 

|cp(6)l  - fl.-'TnTT  u ~f  :r 

m=l  ' mi  i 


(3-53) 


which  for  small  f> , is  approximately 

I / C \ I -l  t 2. 2.  2 \ P 


|cp(6)|  * 1-  (tj  6 w /24)  Z U M .Y , 

\ i mi  / 


if  | f . | = W/2. 


(3-54) 


The  maximum  of  this  under  the  constraints  of  Eq.  (3-39)  occu'rs, 

for  N = 4 and  L = 6,  when 

(j  - U = 1,  1 ^ m ^ 3 

ml  m2 


P = 1 
m3 

H . * 0 
mi 


, 4 £ m ^ 6 


, otherwise. 
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We  then  have 


|cp(6)  * l-5(rr6w)2/8,  N » 4.  (3-55) 

For  6w  = 1/128,  this  represents  an  error  of  about  2.3  x 10  4 which 
is  negligible. 

3 . 3 Choice  of  Values  of  6 in  Parameter  Fitting 

3a.1  x on  the  results  of  the  previous  subsection  we  will  con- 
fine 6 to  the  range 

1/128  £ 6w  * 1/4.  (3-56) 

It  is  also  necessary  to  select  the  grid  of  points  for  6 within 
this  range.  We  have  somewhat  arbitrarily  chosen  to  start  at 
6w  = 1/4  and  successively  halve  the  value  of  6 to  get  the  next 
trial  value.  In  this  manner,  each  tapped  line  can  be  looked  on 
as  a refinement  and  truncation  of  the  previous  one.  We  therefore 
have  the  final  selection  of  values 

6W  = 1/4,  1/8,  1/16,  1/32,  1/64,  and  1/128.  (3-57) 

3 . 4 Range  of  N^. 

Referring  back  to  Eqs.  (3-18)  and  (3-20),  let  us  write 
v^.m  ,n»3)  =o(m1,m2,m3) 

• exp[-j2nN  6wr-3/2-3/N+(m.  +m  +m  )/N]  }. 

— " I **  j 

(3-58) 

For  sufficiently  small  6 we  can  replace  atm^n^m.)  by  its  limiting 
value  (m^ , m2 , m^ ) . Unless  N^6  is  very  small  throughout  the  range 

where  this  approximation  holds  true,  the  values  of  Y will  change 
as  6-0  unless  N^.6  remains  fixed.  We  therefore  can  assume  that 
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N_  " constant/6  as  b-'O.  (3-59) 

j. 

We  will  use  this  property  to  pick  the  range  of  after  the  first 
trial . 

At  the  other  extreme,  we  observe  that  a change  of  nu  by  1 
causes  a change  of 

j2^(NI+  ni)6w/N 

in  the  argument  of  some  exponential  term.  More  accurately,  if 
we  use  an  arbitrary  lattice  of  M^  measurement  points,  the  argu- 
ment of  the  exponential  changes  by 

j (N£  + ni)6w/(M-l), 

which  for  &W  = 1/4  is 

jTT(N][+  n.)/2(M-l)  . 

This  implies  that  if  |n^+  n_jj>  M-l,  then  one  of  the  exponentials 
is  incurring  more  than  90°  phase  change  between  adjoining  lattice 
points.  We  believe  that  such  a situation  corresponds  to  an  under- 
sampled transfer  function  (in  an  engineering  sense)  and  that 
therefore  the  restriction  should  be  assumed  that 

|N  + n. | s M-l, 

or,  since 

1 * n.  * N, 

w-'  will  require  that 

N + 1 * 1-M 

N + N £ M-l 
T 

or 

-M  s N £ M-l-N.  (3-60) 


3-16 


Let  us  now  denote  by  S^(6)  the  apparent  best  choice  of 
for  a particular  value  of  5.  From  Eq. (3-58)  we  would  expect  that 


^(6/2)  « 2*^(6), 


(3-61) 


so  that  the  new  interval  of  delays  covered  by  che  tap  locations 
falls  within  the  old  interval.  We  note  that  with  the  tap  spacing' 

6 the  actual  delays  are  ^ (&)+  1 ]&, . . . , [ft  (6)  + 4] 6 . If  we  then 
choose  the  minimum  and  maximum  values  of  ^(6/2)  to  cover  all  pos- 
sible overlaps  of  this  range  we  must  have 

[Min  N ( 6/2)  + 13(6/2)  - [ft  (6)  l]  6 

"*■  JL 

[Max  Nx  (6/2)  + 43(6/2)  = [^(6)  + 436 

so  that 

2Nj(6)  + 1 £ Nt(6/2)s  2ftI(6)  + 4 .(3-62) 

3 . 5 Best  Choice  of [P  (H)  j 

We  can  now  assume  that  6 and  N^.  are  temporarily  fixed  at 
some  trial  value;  the  values  of  (3 (n)  3 must  be  calculated.  We 
have  available  the  lattice  of  measurements  Y (m  ,m  ,m  ) forlorn.  £ M, 

1 ^ ^ 3.  It  will  be  a labelling  convenience  both  computationally 

and  to  derive  the  best  choice  of  the  { P (n) } if  we  define  two  in- 
dexing variables 


2 

M — M (m1~l)  + M(m2-1)  + m3 

v = N2(n^-1)  + N(n2-1)  + n^. 

The  ranges  of  these  variables  are 

1 £ u £ M3; 

1 £ v s N3. 


(3-63) 

(3-64) 


(3-65) 

(3-66) 


We  modify  the  measurement  values  by  the  first  complex  ex- 
ponential factor  in  Eq. (3-18)  by  setting 


S (u)=Y  (m^,m2,m3) exp| j2TTNI6w|^-  -j  + 


r i+ 

"V  2 "‘3  -1  1\ 

t 2 + 

M-l  J) 

]}• (3-67) 
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We  further  define  a coefficient  array  equal  to  the  triple 
product  of  Eq. (3-18) : 

m.'-l 

C (d,v)  = exp^-jaTT&wVn^-  - + (3- 


We  now  want  to  find  that  se  . of  P (v)  for  which 


V = 


E |S(U) 
U 


E P(v)  C(p,v)  | 

v 


is  minimum. 


(3- 


We  can  interpret  the  array 

C = tc(u,v  ) } 

3 3 

as  an  M -row  by  N -column  matrix;  the  array 

3 = {3  (v)  J 

3 

as  an  N -entry  column  matrix?  and  the  array 
S = (S(U)  ) • 

as  an  M^-entry  column  matrix.  We  can  then  rewrite  Eq. (3-69) 
the  matrix  equation 


where 


V = (S  - Cp)*T(S  - cP)  (3 

* = conjugate 
T = transpose. 


In  expanded  form  this  is 

*T  *T  ♦n*  *P  *T 

v = s s-p  c‘s-scfl  + P c cp.  (3 


Let  us  now  conjecture  that  we  can  rewrite  V as 


*nr> 

v = (p-a)  c C (P-a)  + u 


(3 


where  U is  a scalar  constant  and  a is  a column  matrix  with  w 
entries.  Expanding  the  product  we  have 

*rp  Jf-n  ±rn 

v = 3 c cp-a  c p-p  -c  ca 


*<p  *<p 

+ a c ca  + u. 


(3 


68) 


69) 


as 


-70) 


-71) 


-72) 

3 


-73) 
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If  we  equate  corresponding  terms  in  Eqs.(3-  1)  and  (3-73)  we 
find  that 


so  that 


and  that 


so  that 


*<T»  p 

c ca  = c s 


*T  -1  *T 
a = (C  C)  C S; 


*T  *T  *T 

a c ca  + u = s s 


*T  *T  *T 

u = s s -a  c ca 


(3-74) 


(3-75) 

Substitution  of  the  value  for  a from  Eq. (3-74)  in  this  last  ex- 
pression gives 

* *T  *T  -1  *T  *T  —1 

U = S S - S C(C  C)  C C (C  C)  c s 


or 


*r n 

U = S S - S C(C  C)  c s. 


(3-76) 


This  last  result  can  also  be  phrased  as 

*T  *'T> 

U - S S - S Ca. 


Clearly,  the  minimum  value  of  V is  attained  when 
3 = a 


(3-77a) 


where 


V = U. 


( 3-77b) 


3.6  Program  Outline 

We  are  now  in  a position  to  summarize  the  programs  for  deter- 
mining the  model  parameters. 

3.6.1  Input  Data  from  Measurements 

If  Lhe  input  data  are  to  be  acquired  by  measurements  it  is 
necessary  to  specify 
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a)  the  size  of  the  measurement  lattice  M, 

b)  the  three  carrier  frequencies  V^'V2'  an<*  V3 
(where  is  negative) , 

c)  the  bandwidth  W. 


The  output  of  the  measurements  should  be  a tape  record  or 
card  deck  containing  the  following  data. 

a)  lattice  size  M 

b)  bandwidth  W 

c)  a three  dimensional  array  Yfm^.m^m-j)  where 


Y(m1#m  ,m  ) 


!-W/2  + (m^-1)  W/ (M-l) , 
-W/2  + (m-l)  W/ (M-l) , 
-V7/2  + (m3-l)  W/  (M-l)  J 


(3-78) 


where 


H (v  +f  , v +f  , v +f  ) 

p / f f f \ O -L  -L j 

^l'W  H,  (v^v.+v^+f.+fj 


and 


G = Max  G ( f , f , f ) 

° f f f 12  3 

fl't2'£3 


(3-79a) 


d)  the  normalizing  constant  Gq. 

3.6.2  Input  Data  by  Computations 

If  the  input  data  are  to  be  determined  by  computer  analyses 
of  the  nonlinear  circuit  it  is  necessary  to  duplicate  the  same 
type  of  output  as  in  Section  3.6.1. 

Using  SIGNCAP  for  example,  we  must  specify 
s)  ths  X h 1 1 i C0  cl i rr.0 n s 2. on  M 

where 

2 ^ M £ 5.  (3-80) 
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b)  the  bandwidth  W 

c)  the  carrier  frequencies 

v ,v  .v  where 

v3  < 0. 

One  then  calculates  the  three  dimensional  array  variable 


(3-81) 


Y(m1,m2,m3)  from 


H3l'vl+  flm^),  v2+  f(m2),  v3+  f(m3)] 


' i'  2 * H^+v  +v  + f(m1)  + f(m2)  + f(m3)] 


C3-82a) 


where 


f(mi)  = - W/2  + {m.-l)W/(M-l) 


by  setting 


with 


for  i = 1,2,3;  nu  = 1,2,...,M, 


Y(m1,m2,m3)  = — G(m1,m2,m3) 
o 


G = Max  G(m  ,m  ,m  ) . 
° alllm.)  1 2 3 


(3-82b) 


(3-83a) 


(3-83b; 


Note  that  the  maximum  data  generated  by  this  program  are  125  com- 
plex numbers.  The  program  output  should  then  be 

a)  the  lattice  dimension  M 

b)  the  bandwidth  W 

c)  the  array  variables  Y(m^,ro2,m  ) 

d)  the  normal izer  G . 

o 

3.6.3  Computation  of  Model  Parameters 

The  input  to  the  program  consists  of  the  output  of  the  lat- 
tice value  computation  routine  or  the  measurement  results.  Specif- 
ically the  input  data  should  contain 


3-21 


a)  the  lattice  dimension  M 

b)  the  bandwidth  W 

c)  the  array  variables  Y(m^#m  ,m  ) 

for  1 •'•m.  '•  M;  i = 1,2,  and  3. 

i 

d)  the  normalizer  G . 

o 

In  addition 

e)  the  number  N of  taps/coordinate  must  be  specified. 
This  number  must  satisfy 

1 £ N £ M-l.  (3-84) 


We  define  a normalized  tap  spacing 

D = W&,  (3-85) 

and  take  the  "group  delay"  N^.  as  defined  before.  For  each  D we 

define  the  minimum  value  of  NT  by  N (D) , and  the  maximum  value  of 

I m 

N by  N (D) . We  will  reserve  the  symbol  V for  the  final  minimum 
I M 

value  of  the  approximation  error  and  use  U for  the  approximation 
error  for  a specific  and  D: 

U(K][,D)  = (S-Cp) *T(S-Cp) . (3-86) 

It  will  be  helpful  to  define  partial  minimizations  of  U 


and 


V (N  ,D)  = U(K_,D)  = Min  U(n  ,D),  (3-87) 

T I i n^s  Nj  I 

V (D)  = U (K  , D)  = Min  U(n_,D), 

T o n < Nm  I 

I M 

V_(D)  = U(K_,Dd)  = Min  U(N_,d).  (3-88) 

R R R HSn  I 


Note  that  J<T  is  implicitly  a function  of  N^  and  D,  is  implicitly 
a function  of  D,  and  that  K_,D_  are  implicitly  functions  of  D. 

K K 
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We  can  then  write  the  recursion  relations 


For 


N (D)  = 2K  (2D)  + 4, 
M o 

N (D)  = 2K  (2D)  + 1; 
m o 


( 3-89a) 
( 3-39b) 


VVD) 


j Vt(Ni-1,D)  if  U(NI,D)  > V(p(NI-l,D)^ 

V U(N;[.,D)  if  U(NI#D)  ^ Vt(N];-1,D) j 

(3-90) 


K (N 
T I 


Km(Nx-l,D)  if  UtNj.D)  > Vt(Ni-1,D) 


,D)  . f VV— 

^ N_  if  U (N_ , 


rD)  s Vt(Nj-1,D) 


vt(d)  = vt(nm,d), 

K = K (N  , D)  ; 
o T M 


)■ 


(3-91) 

(3-92) 

(3-93) 


Vr(D)  = 


Vr(2D)  if  Vt(D)  > Vr(2D) 
Vt(2D)  if  Vt(D)  * vr(2D) 


(3-94) 


kr(D) 


■{ 


KR  (2D)  if  V,p(D)  > Vr(2D) 

K (d)  if  V (D)  * V (2D) 
O T R 


}• 


(3-95) 


VD)  - 


DR (2D)  if  Vt(D)  > VR (2D) 


D if  V (D)  * Vr(2D) 


(3-96) 


For  each  of  these  recursions  we  also  need  initial  conditions, 

N and  N we  have 
m M 


{ N = - M 

1 N = M-l-N 
V.  M J 


when  D = 1/2, 


(3-97) 
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n ri77^’^T77i,' 


Since  we  have  constrained  the  values  of  y by 

2 

|y(m  ,m  ) | ^ 1, 

it  follows  that 
*T 

S S £ 125 

so  that  V and  V can  be  initialized  by  setting  them  equal  to  any 
R T 

number  larger  than  125,  say  200.  We  have  then 

v (N  -1,D)  = 200  (3-98) 

T m 

V (1/2)  = 200. 

R 

The  final  values  of  the  parameters  then  yield  the  optimum 
parameter  values: 

Optimum  value  of  6w  = D (1/128)  (3-99a) 

Optimum  value  of  N = K (1/128) 

1 R 

Minimum  of  V = V (1/128) . 

K 

We  then  have  the  following  scheme  for  calculations: 


V/Siy 


Note  that  under  the  program  restrictions, 
1 s U i M3  s 125. 


e)  For  each  nu, 

1 £ m.  ^ M;  i = 1,2,3 

and  for  each  n. 

i 

1 ^ n £ N;  i = 1,2,3, 
define  M as  before  and  define 

v = Nirij-  1)  + N (n2~  1)  + n^, 

and  a two  dimensional  array  variable 


c(m3v) 


1 

2 


+ 


m.  - 1 
i 

M-l 


Note  that 

3 

1 * u s M * 125, 

1 s v £ N3  s (M-l) 3s  64. 

f)  Interpreting  S(u)  as  a column  matrix  and  C(n,v) 

3 3 

as  a rectangular  matrix  with  M rows  and  N 
columns,  calculate  the  number 


*T  -1  *T  , *T>  *T  -1  *T 
U=Cc(C  C)  C S-S]  [C(C  C)  C S-Sj 

g)  If  U < V , 


Make  available  for  print  out  N^ , D,  and  U, 


h) 

Replace  NT  by 

JL 

Nt+  i 

X 

i) 

If  N * N go 
1 M 

back  to  step  (d) 

j) 

Set  K = K 
o T 

T* 


j 
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k) 

l) 


m) 

n) 

o) 


P) 


q) 


r) 


s) 


t) 


If  Vm  < V set  {V=  V , K = K , D = d}. 
T R R T R T R 


Set  N = 2K  + 5 
M o 


N = 2K  + 1 
I o 


D = D/2 


Vm  = 200. 

T 


If  129  D ^ 1 go  back  to  step  (d) 


Set  D = Dr,  N][  = Kr.  V = V . K - Kr 


Calculate  the  array  variables  S and  C as  in 
steps  (d)  and  (e) 


Calculate  the  array  variables 

3 

3(v)  for  v = 1 to  N by  interpreting  0 as  a 
column  matrix  given  by 


*T  —1  *T 

0 = (C  C)  C S. 


for  1 ^ n^  £ N,  for  i = 1 to  3 set 


B*nl'n2'n3*  =^[(i11"1)n2+  (n2~  1)  N + 
Calculate  the  normalized  error 


P = 


V 


*T 

S S 


Make  available  for  print  out 


GQ, r D, p and  the  array  variables  B(n^,n 


Make  av  able  on  cards 


Gq/K,6  = D/W,N,  and  the  array  variables 


B (n^n  ,n  ) . 


^^AI%<^JCjjrITCTV  ^(Ktv  *■ 


SECTION  4 

SAMPLED  DATA  BUTTERWORTH  FILTER 

4.1  Use  of  Approximate  Linear  Filter  Response 

After  adding  the  third  order  distortion  terms  to  the  linearly 
amplified  components,  it  is  necessary  to  filter  the  resultant  sig- 
nal so  that  its  bandwidth  is  equal  to  the  nominal  IF  bandwidth  of 
the  receiver.  Simulation  of  the  exact  IF  response  on  the  computer 
would,  in  general,  require  a convolution  involving  a very  large 
number  of  delayed  replicas  of  the  input  process,  and  would  be  very 
time  consuming. 

It  will  be  adequate  for  most  purposes  to  replace  the  impulse 
response  of  the  actual  linear  circuit  by  any  filter  having  the  cor- 
rect 3 dB  bandwidth  and  adequately  fast  roll  off.  It  should  be  em- 
phasized that  such  a change  does  not  have  the  drastic  effect  on 
nonlinear  spectra  that  would  result  from  changing  the  actual  IF 
transfer  function.  The  effect  of  this  actual  receiver  filter  will 
have  been  incorporated  in  the  equivalent  transfer  function  K(f^,  f^)  ; 

the  substitution  of  an  approximate  filter  for  H(f)  causes  only 
linear  distortion  of  the  output. 

The  filter  that  is  chosen  to  approximate  H(f)  should  obviously 
be  selected  for  ease  of  simulation  on  the  computer  as  well  as  for 
reasonable  match  to  IF  filter  characteristics.  The  second-order 
Butt_^worth  filter  meets  all  these  requirements  satisfactorily:  it 
has  maximally  flat  inband  response,  it  yields  12  dB  rejection  at 
one  bandwidth  separation  from  its  center  frequency,  and  it  can  be 
simulated  by  a two  stage  recursive  filter  on  the  computer. 
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4.2  Use  of  Butterworth  Filter  for  Spectral  Shaping  of  Noiselike 
Signals 

In  order  to  simulate  the  effects  of  random  interference  it 
is  desirable  to  include  a noiselike  waveform  as  one  of  the  possible 
inputs.  This  waveform  can  be  generated  conveniently  using  a sequent 
of  random  complex  numbers  at  the  sampling  rate  of  the  simulation. 
However,  it  is  necessary  Lo  provide  spectral  shaping  in  order  to 
approximate  the  bandwidth  characteristics  of  the  interference. 

We  have  again  chosen  a second-order  Butterworth  sampled  data  filter 
for  chis  function. 

4.3  Properties  of  Second-Order  Filter 

Let  {z^:  £ 00 } be  a sequence  of  complex  numbers  derived 

from  sampling  an  input  process  using  a sampling  interval  6,  and 
let 

t = (2PcosCp)  Y - p2  Y _+  AZ  , (4-1) 

n n-i  n-2  n 

where  P and  cp  and  A are  real  constants.  We  assume  the  nonrecursive 

definition  of  Y to  be 
n 


Y = XL  (a.b  m + a b_m)  Zn  t 
n m— 0 11  2 2 n-i 

Then  since  Eq.  (4-1)  can  b written  as 


(4-2) 


Y - (2pcosCp)  Y .+  p Y _ = AZ 
n n-1  n-2  n 


(4-3) 


we  can  determine  a^*b^,a2,b  by  substitution  of  Eq. (4-2)  in  Eq. 
(4-3)  to  yield 

CD 

<V  a2>V  (a!bl+  a2b2)Z„-l+  Kio(albl+k+  a2b2+k)Zn-2-k 

00 

- (2PCOSCP)'  (a  + a ) Z .+  Z n(a.b.i+k+  aj>^1+k)z  0 ,1 

L 1 2 n-1  k=0\  11  22)  n-2-kJ 

+ p2|"  E (a,b  k+  a b k)z  0 1 

Lk=0  \ jl  1 2 2/  n-2-kj 


)z  1 
/ n-2-kj 


= AZ 
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fc-ts ,»>«#=>  ^wjkks*?*???**’  y,  *sj5*s f^s "nsv  y <v<r- ;-  - ;*- ■ <v  * '«-* 


From  this  equation  we  deduce  that 

2 2 
b.  - (2pcoscp)b.  + P = 0,  i = 1,2; 

l i 


al  + a2  - A; 


albl+  a2b2~  (2pcosCp^  ^ai+  a2^  = °* 


(4-4) 

(4-5) 

(4-6) 


Equation  (4-4)  has  the  solution 


r 2 2 . 2 

2pcoscp  a-  y4P  cos  cp  - 4P 
bl  = 2 


or 

b^  = P exp  (jcp) 
b2  = P exp  (-jcp)  . 

We  observe  that  a2  — a^  since  b2=  b]_  * can  thus  write 

a^  ~ I a | exp  ( ja)  , 
a2  = | a | exp  (~ja)  , 
so  that  Eqs . (4-5)  and  (4-6)  become 
2 | a | cos  a = A 
— 2 1 a | P cos(P-a)  = 0. 

We  deduce  immediately  that 
a = cp  - rJ/2 


and 


Hence 


|a|  =■  A/(2sincp)  . 


a^  = (a/2)  (1-jcotcp)  , 


/ Jk  /O  \ { 1 -i--i 

\ r\/  a-  ; 


(4 -7  a) 
(4 -7  b) 

(4 -8  a). 
(4 -8b) 

(4-9) 

(4-10) 


(4 -11a) 


IA- 1 


We  thus  have 


Y = (A/2)  1 0m  Z r(1"jCOtCp)  exp  (jniCP)  1 . (4-12) 

Yn  m=0  n-m  L+(l+jcotcp)  exp  (-jmcp)J 
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4.3.1  Response  to  Sampled  Cissold 

Now  let  [z  } be  the  samples  of  a complex  cissoid  with  fre 
m 

quency  f: 

Z = exp  (j2^f6m)  . (4 

m 

We  can  shorten  subsequent  expressions  by  setting 


P = 2"f6  (4 

so  that 

Z - exp  (jmP)  . (4 

Substitution  of  this  in  Eq.  (4-12)  yields 

2Y^  exp(-jn(3)/A 
00 

= (l-jcotcp)mEQ  [p  exp  (jcp-jp)]m 

GO 

+ (l-jcotcp)in|:0  [p  exp  (-jcp-jP)] 


or 


so  that 


1-jcotcp  l+jcotcp 

l-pexp(  jcp-j3)  l-pexp(-jcp-jP)  ' 

2 jY^  (sincp)  exp  (- jnp  ) /A 

= [exp (- jcp)  - p exp  (-jP)]  1 
-Texp(jcp)  - p exp(-jP)]  1. 

2 -1 

=2jsincp  ri-(2pcoscp)  exp(-jP)  + P exp(-j2P)] 

2 -1 

Yn  = Aexp(jnP)  [1- (2pcoscp)  exp(-jP)  + P exp(-j2P)] 

(4 


We  then  have  the  magnitude  of  Y given  by  the  expression 


| | 2=  A2  1 1-  (2pcosfp)  exp  (- jP)  + p2  exp(-j2p)j 

= A2  | [exp(jcp)-  p exp(-jP)  ] [exp(-jcp)-  p exp(-jP)]j  2 
= A2  f l-2pcos  (P+cp)  + P2]  1[l-2Pcos  (P-cp)  + p2]  * 

| Y |2/A24"U+  P2)2  - 4p  ( 1+P 2 ) cos  P cos  t P + 2p^  (cos  2P+cos2cp)~| 


(4-17) 

If  we  now  write 

cosP  - 1-2  sin2y  (4-18a) 

where 

Y - P/2,  (4-18b) 

and  then  use  the  relation 

2 2 4 

cos  28  - 2cos  8-1  = 1-8  sin  y + 8 sin  Y,  (4-18c) 

Eq.  (4-17)  can  be  rewritten  as 

A /Iy  I'  = (1+p  ) 4P  (l+P^)  coscp  + 2P  (l+cos2cp) 
n 

2 2 2 2 2.4 
+8p  (1+P  ) coscp  sin  Y - 16p  sin  Y + 16p  sin  Y 


or  as 

A2/Iy  I2  = (l-2pcos  c?+p^)2 
n 

2 2 
+ 8p  f (l+p  ) coscp-2p]  sin  Y 

2 4 

+ 16p  sin  Y.  (4“19) 

4.3.2  Special  Choice  of  P 

2 

We  now  choose  P so  that  the  coefficient  of  sin  Y vanishes 
in  order  to  have  a maximally  flat  response: 

2 

(1+P  ) coscp  - 2P  = 0 
or 

2 

cos  cp  = 2p/  (1+P  ) . 


(4-20) 


We  then  have 


sincp  = (l-p2)/(l+p2) 


(4-21) 


so  that 


2 _ 1-sincp  _ 1-cos  (cp-rr/2) 
1+sincp  1+cos  (cp-tt/2 ) 


(4 -22c 


or 


= tan  (co/2  - tt/4) 


p = tan  (cp/2  - ”/4)  . 


(4 -22b 


With  this  choice,  the  leading  term  of  Eq.  (4-19)  becomes 

2 2 

(1-2P  coscp  + p ) 


-Ti- 


4p 


L 1+p2 


+ P 


21 


o 2 2 

= r(izO  i 
"L  i+p2- 


Hence 


or 


. 2 2 
~2sin  cp~] 

|_l+sincpj 


. 2 J2. 


i.2 1 ,L  .2  T2sin  cpT.  . ,ri-sincpl  . . 

I A / hr/  = T7 — : — 1 + 16  — — 'sxnY 
1 n1  Ll+smcp_,  Ll+smcp_j 


2 4 -1 

|y  |2  2ri+sinco  *j  p.  4cos  cpsin  Y~1 

n L . 2 _ u .4  J 

2sin  cp  sm  cp 

In  order  to  have  unity  gain  at  zero  frequency  we  take 


(4-23) 


A = 


„ • 2 
2sin  cp 

1+sincp. 


(4-24) 
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Furthermore,  the  3 dB  attenuation  point  is  achieved  when 

2 ' 4 

4 cos  cp  sin  Y 

.4  " 1 (4-2! 

sm  cp 

If  we  wish  this  3 dB  point  to  be  at  f = w/2,  then  when  0 = nw6 
(see  Eq.(4-14) ) or  when  Y = nW&/2  (see  Eq.  (4-18) ), Eq. (4-25)  must 
be  satisfied.  We  thus  require 

• 4 

sm  cp  . 4 

2 = 4 sln  Y0  (4-2C 

COS  cp 


where 


Y = ttW6/2  . 
o 


We  can  solve  Eq. (4-26a)  immediately  to  arrive  at 


cp  = arc  sin  s 


where 


s-r 


r 2sin2Vo  1/2 

Lri+sin4Y  ]1/2  +°sin2Y  ->  ' 


We  thus  have,  upon  substitution  in  Eq.  (4-24) , 
A = 2S2/(1 +S), 

and,  upon  substitution  in  Eq.  (4-22a), 

P = C(l-S)/(1+S)]1/2. 

Finally,  since 

cos  cp  = (l-s2)1//2  = [ (1+S)  (1-S) 
the  recursion  relation  Eq.  (4-1)  can  be  written  as 


r / 1 / / -1  . \ 1 


= 9 / 1 _C\  \r  r / 1 n\  / /1  \ . r 2 , 

An-1  ~ \a-o//  j x ti  zt>  / U+SJ  JZ  , 

11  -1-  n—  ^ n 


(4-31) 
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4.3.3  Response  to  Independent  Samples 


We  now  investigate  the  response  of  this  sampled  data  filter 
to  a sequence  of  independent  complex  samples.  Repeating  Eg . (4-12) 
here, 

Y = (A/2)  E om7  rU-jcotcp)  exp  ( jmcp)  1 U-32) 

n m=0  n-m  L+(l+jcotcp)  exp  (-jmcp)J 


Letting 

Vz  = E{|znl2} 

be  the  common  variance  of  the  input  sequence,  we  have 

2-;  . 2. 


VY  = E(|Yn!  } " (A  V4) 


(4-33) 


*F.  P 

m=0 


2m 


(1-jcotcp)  exp  (jmcp)  i2 
+ (l+jcotcp)  exp  (-jmcp) 


(4-34) 


We  can  immediately  rewrite  this  as 
4 (sin^cp)  Vy/A^V^ 

CO 

= m£0  P I exp  ( jcp+jmcp)  -exp(-jcp-jmcp)  I 


= - exp  ( j2cp)m|of(?exp  ( j2cp)  ] 


m 


00  2 

- exp  (- j2cp)  Cpexp  ( — j 2cp)  ] m 


+ 2 S p 
m=0 


2m 


exp  (-j2cp)  __2_ 


= 2 


exp  (j2co) 

2 **  2 ^ 2 
1-p  exp  ( j 2cp ) 1-p  exp  (— j 2cp ) 1-p 

(1+P)  [1-cos  (2cp)  ] 


2\r  1 JL  a 

\ -t  “ K / L J-  * K 


4 2 /0  , , 

O/s  nnc*  I 7#A  I 1 

— ah  / -j 
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If  we  now  substitute  Eqs.  (4-27)  through  (4-30)  in  this  we 


obtain 


V. 


V. 


(4-35) 


In  order  to  provide  a unit  variance  output,  we  therefore  set 


V_ 


2-S 


(4-36) 
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SECTION  5 


NONLINEAR  TRANSFER  FUNCTION  FOR  TESTING  COMPUTER  ROUTINES 


It  is  helpful  to  have  available  a routine  for  generating  the 
first  and  third  order  r.ransfer  functions  H^(f)  and 
that  does  not  require  the  use  of  SIGNCAP.  This  section  presents 
an  outline  of  a routine  for  providing  such  transfer  functions. 


Consider  the  circuit  diagram  shown  in  Fig.  5.1,  consisting 
of  two  single  tuned  circuits  separated  by  an  amplifier  that  ex- 
hibits a cubic  distortion  term.  If  H (f)  is  the  transfer  function 

A 

of  the  first  filter  then  an  input  of  the  form 

j2rrvnt 


v (t)  = E a e 
o n n 


(5-1) 


where  the  occur  in  pairs  of  positive  and  negative  frequencies, 

yields  an  input  to  the  nonlinearity  of  the  form 


v (t)  = >1  a H (v  ) e 
1 n n A n 


j2nv  t 
n 


(5-2) 


The  output  of  this  nonlinearity  is 


- ; a H (v  ) e 
5 n n A n 


j2nv  t 
n 


+ B l 


ni  n2  n3  ni 


i.  a H (v  ) H (v  ) H (v  ) 
n2  n3  A n;L  A n2  A n3 


exp 


r j2rrt  ( 


V + V + V ) 

ni  n2  n3 


(5-3) 


The  second  linear  filter  has  the  same  transfer  function  as  the 
first  no  that  the  final  output  is 
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2 


v (t)  = S a H (v  ) e 
3 n n A n 


j2nv  t 


n 


+ B T T.  >:  a a a H (v  ) H (v  ) H (v  ) 
ni  n2  n3  "l  n2  n3  A ni  A n2  A n3 


H,  (v  + v + v ) 

A n2  n3 

exp^j2nt(vn  + + vn  ) i.  (5-4) 


Since  we  can  also  write 


v. 


(t)  = J< 


df  e 


j2nft 


H.  (f)  E a 6 ( f-v  ) 
Inn  n 


j2TTt(f1+  f2+  f3) 


'2  3 


WW 


E E >:  n a 6 (f . -v . ) 

ni  n2  n3  n-  1 i 


(5-5) 


i=i 


we  have  the  immediate  correspondence 


Vf)  - H^(f) 


VWV  = BHA(V  f2+  f3}  fi  W* 

i=l 


(5-6) 

(5-7) 


The  transfer  function  of  the  single-tuned  filter  is  just 

«A(f)  = [1  + jR[2nfC  - 1/ (2TTfL)  j }_1 . (5-8) 


and 


If  we  use  the  usual  notation 
fQ  = l/2n (LC) 1/2 , 


U = K/  h,‘  , 

o 


(5-9) 


(5-10) 
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I 


we  can  write  the  transfer  function  as 

-1 

H (f)  * [1  + jQ(f/fQ  - f/f)]  * (5~n) 

For  numerical  calculations  to  test  various  program  routines 
it  is  necessary  to  assign  numerical  values  to  Q and  fQ,  and  also 
to  select  the  bandwidths  and  center  frequencies  of  the  signals  to 
be  accommodated  by  the  canonic  model.  To  this  end  we  take 


f = 50  MHz 
o 

°r  7 

fc  = 5 x 10 

and 

Q = 102 


(5-12a 

(5-12b 


to  characterize  the  filter.  We  will  then  model  a crossmodulation 
situation  where  the  desired  carrier  frequency  v^  is  at  the  center 
frequency  of  the  filter: 

v = f = 5 x 107,  (5-13a 

1 o 

and  where  the  interfering  carrier  is  1 MHz  removed  from  the  desired 
carrier: 

v2  = 5.1  x lo7  (5-13b 

v = - 5.1  x 107.  (5-13c 

We  will  take  the  nominal  bandwidth  of  the  signals  to  be  0.5  MHz: 


W = 0.5  x 10 


6 


(5-13d 


These  choices  of  reasonably  realistic  transfer  functions  and 
numerical  values  will  permit  testing  of  the  parameter-fitting 
routines  without  the  necessity  of  time  consuming  calculations  with 
nonlinear  circuit  analysis  programs. 


« 

It  is  an  additional  convenience  in  testing  program  routines 
to  include  the  effect  of  computational  or  experimental  error.  This 
can  be  accomplished  by  introducing  deliberate  round-off  error  in 
the  or  outputs. 
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SECTION  6 


COMPUTER  GENERATED  WAVEFORMS 

After  having  computed  the  model  coefficients  by  the  procedure 
outlined  in  Section  3,  the  necessary  parameters  are  available  for 
calculating  the  complex  envelope  of  a third  order  interaction  term. 
The  general  form  of  this  complex  envelope  can  be  found  by  substitu- 
tion of  the  tapp-  i delay  line  model  into  Eq. (2-7)  to  yield 


a (t)  2 f3(n)  ft  z (t  - n.6  - K6) 
p n — . 1 1 n . i 

& — i=l  Fi 


(6-1) 


where  K is  the  optimum  value  of  found  by  the  procedures  out- 
lined in  Section  3.  For  purposes  of  computation  we  will  actually 

evaluate  a (t)  only  at  integer  multiples  of  6,  however.  Let  us 
£ 

therefore  introduce  the  notation 


and 


Y_(J)  = a (J6  + K6) 

3 R 


Z. (J)  =z  ( J6) . 

1 pi 


(6-2) 


(6-3) 


(It  will  be  noted  that  these  represent  a duplication  of  earlier 
use  of  upper  case  letters  for  spectra;  no  confusion  should  result 
in  the  balance  of  this  material  where  only  sampled  data  sequences 
are  to  be  considered.)  We  then  have 

3 

Y3(J)  = 2 p(N)  jf[  Zi(d  - N ).  (6-4) 

N i=l 

This  sampled  output  Y^fJ)  is  advanced  by  K samples  with  respect 
to  the  linear  output  term  and  this  discrepancy  of  alignment  must 
be  compensated  in  the  program  by  delaying  Y^  by  K samples  before 
combining  if  K is  positive,  or  delaying  the  samples  of  the  linear 
output  by  I K | samples  if  K is  negative. 
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It  will  also  be  recalled  that  in  Interference  Cases  #1  and 
#4  of  Section  2.6,  {Eqs. (2-34)  and  (2-37)),  the  interference  en- 
velope is  actually  modulating  the  adjacent  carrier  frequency  sepa- 
rated by  W Hz  from  the  linear  output.  Hence  a transformation  of 
the  samples  tY^(J)J  equivalent  to  this  frequency  translation  must 
must  be  affected. 


Finally,  as  discussed  in  Section  2.7,  the  amplitude  of  the 
third  order  product  must  be  scaled  relative  to  the  linear  component 
to  take  into  account  both  the  relative  amplitudes  of  interferers 
and  desii ed  signal,  and  to  incorporate  the  normalizing  constant  Gq 
of  Section  3. 


Let  {Y(J)}  be  the  samples  {Y^ ( J) } corrected  for  the  frequency 


offset : 


Y ( J)  = e2,rAW(;J  2 B(N)  TT  Z.  (J-N.  ) 
N _ i 1 1 


(6-5a) 


where 


fO  if  interference  product  is  at  nominal  carrier  ") 

A ={  } 

1^1  if  interference  product  is  at  adjacent  carrier  J 

(6- 5b) 


and  let 

Y„  (J)  = Z (J+K)  (6-6) 

u o 

with  {z  (J) ] the  samples  from  the  desired  signal.  We  can  then 
write  the  samples  of  the  total  output  from  the  unit  gain  amplifier 
and  idealizied  third  order  transfer  function  as 


Yt(J)  » ClYl(J)  + CtY  (J)  (6-7) 

where  C and  C are  normalizing  constants  that  include  the  con- 

Xj  j; 

scant  and.  the  relative  powers  of  the  signal  and  interferers 
as  well  as  the  overall  gain  of  the  receiver.  For  most  computations 


* 


where  only  relative  distortion  is  important,  only  the  ratio  of 

C to  C need  be  specified,  and  the  absolute  scaling  can  be  ac- 
T L 

complished  for  computational  convenience. 

6 . 1 The  Four  Siqnal/Interferonce  Combinations 

In  Section  2 we  discussed  the  four  cases  appropriate  for  anal- 
ysis. We  now  list  these  cases  in  the  notation  of  this  section,  using 

the  additional  notation 

S(t)  = signal  waveform  -x 

U^(t)  = interfering  waveform  \ • (6r8) 

U (t)  = different  interfering  waveform  J 
2 

Table  6.1  Signal/Interference  Combinations 


Case  1 

Case  2 

Case  3 

Case  4 

z (t) 
o 

S(t) 

S(t) 

5 (t } 

S(t) 

z ^ (t) 

ux(t) 

ux  ( t) 

S(t) 

ux(t) 

z2(t) 

u1(t) 

ux(t) 

ux(t) 

U2(t) 

z3  (t) 

ux*(t) 

u2(t) 

U *(t) 

V<t) 

A 

1 

0 

0 

i 

Eauivalent  IF  Filter 

It  is  necessary  to  include  the  effect  of  the  linear  filter 
H(f)  of  the  equivalent  receiver  to  determine  the  overall  impact 
of  the  interference  on  inband  interference.  As  we  discussed 
earlier,  it  is  adequate  to  approximate  this  filter  by  any  filter 
which  restricts  the  bandwidth  to  W Hz;  and  a sampled  data  second- 
order  Butterwortn  filter  is  adequate  for  this  purpose. 
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Using  the  results  of  Section  4,  we  define 


y = ttw6/2 


(6-9) 


S = ! 


2 sin  y 


■1 


,1/2 


L,,  .4,1/2  . 2 J 

(1+sin  Y)  + sxn  Y 

The  filtered  output  sequence  is  then  given  by 

2 

1_C  OC 

Y (J)  = 2 (1-S)  Ymri  (J-l)  - TTT7  Y ( J-2)  + ~ Ym(J). 
TF  TF  1+S  TF  1+S  T 


(6-10) 


(6-11) 


We  can  also  examine  the  filtered  version  of  the  third  order  produc 
without  the  linear  term  by  defining 

Yp(J)  = Ytf(J)  when  CL  = 0.  (6-12a 

Correspondingly,  we  can  define  a linearly  filtered  signal  by 

YLF(J)  = Ytf(J)  When  CT  = °*  (6"l2b 


6. 3 Signal  and  Interference  Waveforms 

The  model,  as  it  stands,  will  accept  any  choice  of  waveforms 
for  the  signal  and  inter ferers  of  Table  6.1.  It  is  a program  re- 
quirement that  it  be  possible  to  generate  typical  waveforms  in- 
ternally during  the  computations.  We  believe  that  suitable  wave- 
forms can  be  provided  by  the  following  repertoire: 

a)  Signal  waveform: 

i)  CW 

ii)  Sinusoidal  amplitude  modulation 

iii)  Sinusoidal  frequency  modulation 

iv)  Combined  FM/AM  with  different  modulating  frequencies 

b)  First  Inteferer: 

Same  possible  characteristics  as  for  signal 
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c)  Second  Interferer: 


Noiselike  waveform  with  second-order  Butterworth 
frequency  characteristics. 

We  therefore  define  the  following  sampled  data  sequences 

ST  (J)  = [1  + M cos  (2tt  f J6  + 9 )] 

IG  S S - w s 

AM 


exp[jD  cos  (2^  J6)]; 

s sml 

FM 


(6-13) 


U(J)  = [1  + H cos  (2^  f J6  + 0)3 

XAM 


expTj  D cos (2™  f J6)]; 

FM 


(6-14) 


2 c S 

n(J)  = 2(i-s)n(J-i)  h (J-2)  + - ~r-  g(j). 


where 


(6-15a) 


= r [1-(1-2S2)  (l-s2)  1/23~] 1/2 
5 ^2s3d+$)ri+(i-s2)1/2] 


(6-15b) 


LG(J)}  are  independent,  zero  mean  complex  Gaussian  random 
variables,  with  unit  variance  (6-15c) 


and  the  initial  conditions  are 


rj(l)  = r,(2)  = 0.  (6-15d) 

(The  normalizing  constant  C_  was  derived  in  Section  3.) 

5 

In  choosing  the  modulating  frequencies  and  deviation  ratios 


it  is  necessary  to  insure  that  the  resultant  bandwidth  of  the  complex 
envelopes  does  not  exceed  W.  We  note  that  a deviation  ratio  of  ap- 
proximately 2.405  permits  generating  an  FM  waveform  having  complete 
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carrier  suppression;  we  therefore  have  chosen  to  restrict  the  maxi- 
mum deviation  ratio  to  2.5: 


0 * DS 
0 £ DT 


* 2*5  ]. 
£ 2.5  J 


(6-16) 


At  this  peak  deviation  ratio,  the  sideband  power  distribution  is  as 
given  in  Table  6.2. 

Table  6.2 

Sideband  Power  for  Deviation  Ratio  of  2.5 


Component 

Carrier 

Sideband  # + 1 

Sideband  # + 2 

Sideband  # + 3 

Sideband  # + 3 


Fraction  of  Total  Power 
« 0 

0.25 

0.20 

0.03 

0.005 


Restriction  of  the  FM  modulating  frequency  so  that  the  fourth 
sidebands  are  included  in  (-W/2,W/2)  is  adequate  to  meet  the  band- 
width restriction.  In  the  AM  and  combined  AM/FM  cases,  the  ampli- 
tude modulation  introduces  an  additional  spreading  of  every  com- 
ponent equal  to  the  amplitude  modulating  frequency.  We  therefore 
need  to  restrict  the  pairs  of  modulating  frequencies  by  some  re- 
lation of  the  form 

f + 4f  <:  w/2  , (6-16a 

SAM  SFM 

f + 4f  ^ W/2.  (6-16b 

*AM  FM 

6.3.1  Drift  Frequency 

Because  of  the  choice  of  sample  i.  *acing  to  be  an  integer  sub- 
multiple of  the  reciprocal  bandwidth,  and  because  of  the  placement 
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of  an  interfering  carrier  exactly  one  bandwidth  away  from  the 
desired  signal  in  Cases  1 and  4 of  Table  6.1,  it  is  possible  for 
a CW  or  AM  interferer  to  yield  an  interference  product  in  a fixed 
phase  relation  relative  to  the  desired  signal.  To  avoid  this  pro- 
gram artifact  it  is  useful  to  introduce  a "drift"  frequency  in 
Eq. (6-14) . This  drift  frequency  should  be  chosen  to  be  less  than 
a few  percent  of  the  bandwidth, and  irrationally  related  to  it. 

These  requirements  can  be?  met  by  taking  the  drift  frequency  to  be 
W/20tt  so  that  Eq.  (6-14)  can  be  replaced  by 

U(J)  = [1  + u cos(2nf  j6  + 0 )] 

*AM 

exp  |j^0.lJ6w  + cos(2TTfi  j6)Jj.  (6-17) 

FM 

6.3.2  Interference  Cases  with  Internal  Routines 

In  Table  6.1  we  listed  the  possible  signal/interference  com- 
binations. With  the  internal  routines  available  for  generation  of 
"modulated"  interference  (the  sequence  U(J))  and  noiselike  inter- 
ference (the  sequence  t(j)),  it  is  possible  to  create  two  distinct 
interference  products  for  each  case,  depending  on  how  we  associate 
and  with  U and  ri.  We  can  therefore  expand  the  table  to  yield 
the  eight  cases  shown  in  Table  6.3. 

6.4  Demodulated  Outputs 

It  is  of  some  help  in  evaluating  the  impact  of  nonlinear  dis- 
tortion to  determine  the  distortion  voltage  after  detection.  It  is 
possible  to  provide  this  option  in  the  computer  routines  for  both 
envelope  detection  and  phase  detection.  It  should  be  observed 
that  for  both  of  these  types  of  detection  the  "filtered"  output 

version  [y  ( J) } or  {Y  (J)  } in  Eqs. (6-lla)  and  (6-12b)  should  be 
TF  LF 

used,  since  the  unfiltered  distortion  products  have  a bandwidth  of 
3W. 
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Table  6.3 


Internally  Generated  Signal/Interference  Combinations 


Offset 

Parameter 


Case 

A 

VJ) 

la 

1 

SIG(J) 

lb 

1 

SIG(J) 

2a 

0 

SIG(J> 

2b 

0 

SIG(J) 

3a 

0 

SIG(J) 

3b 

0 

SIG(J) 

4a 

1 

SIG(J) 

4b 

1 

SIG(J) 

z1(J) 

z2(J) 

z3(J) 

U(J) 

U(J) 

U*(J) 

T|(J) 

n(J) 

n*(J) 

U(J) 

U(J) 

h(J) 

T1(J) 

ri(J) 

U(J) 

SIG(J) 

U(J) 

U*(J) 

SIG(J) 

h(J) 

h*(J) 

U(J) 

h(J) 

•n*  ( j) 

h(J) 

U(J) 

U*(J) 

A simplified  presentation  of  the  distortion  in  the  detector 
outputs  is  possible  if  only  the  distortion  is  made  available  at 
the  output.  In  the  case  of  the  phase  detector,  this  also  simplifie 
the  computation  in  that  it  eliminates  computational  errors  of  2tt 
in  computing  the  arc  tangent.  We  therefore  define  an  envelope  dis- 
tortion 


env(J)  = 1 ytf ( J) I - I ylf ( J ) I 


(6-18) 


and  a phase  distortion 


ph(J)  = arg[YTF(J)]  - arg|YLp(J)|. 


it  Wlii  De  noted  tnau  uua  Acn-eei"  expression  can  we  wi" i e i_en  as 
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ph  = Im  { log(YTF/YLp)} 


= I,  {log  [l  + 

L LF  -J 

which,  if  the  distortion  is  small  enough  to  be  tolerated  at  all, 
can  be  approximated  by 


ph  « Im  j 


Y - Y 
TF  LF] 


LF 


■ Im  • 


We  will  use  this  approximate  formula  as  sufficiently  precise  for 
the  purpose  of  estimating  phase  distortion.  (The  phase  distor- 
tion can  be  equated  to  rr  radians  when  Y is  small.) 

LF 

6. 5 Necessary  Computer  Routines 

6.5.1  Generation  of  Signal,  Interference  and  Noise 

(a)  The  necessary  input  parameters  are  given  in  Table  6.4. 

Table  6.4 

Input  Parameters  for  Waveform  Generation 


Description 

Textual  Notation 

Bandwidth 

w 

Tap  spacing 

6 

AM  modulation  index, signal 

U 

s 

AM  modulation  index, interferer 

FM  deviation  ratio, signal 

Ds 

FM  deviation  Ratio, interferer 

DI 

AM  modulating  frequency , signal 

AM  modulating  frequency, interferer 

ru  i 

flAM 

AM  phase, signal 

P 

S 

AM  phase,  interferer 

6I 

FM  modulating  frequency , signal 

SFM 

FM  modulating  frequency , interferer 

fr 

XFM 
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(b)  Restrictions  on  parameters 

W6  5 i/4 ; 

0 5 V 1/ 

0 * Ui«  1; 

0 * Dg::  2.5, 

0 * D s 2.5 


If 


U = 0 


r 


D.  > 0 


i then 


fS  0 
AM 

0 s f £ W/8 
SFM 


ns  = ° 


J 


If 


Ds  = 0 


US  > ° 


then 


( f S 0 

1 FM 


0 * f * W/2 
AM 


. 


L“s  - o 


J 


If 


( 0 s f 5 W/10  >. 
i AW  I 


o 5 f ^ W/10 
bFM 


0 s P s 2n 

V,  o 
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If 


fDl  - 0, 

lWT  > Oj 


then 


( f =0 

FM 

0 s f <:  w/2  \ 

iAM 

\ = ° 


ro  i f.  £ W/10 

XAM 

It(Dl  > ° "I  then  | s * W/1° 

Ur  -'  0 j 


l 


>• 


0 £ A^.  £ 2" 


(c)  Recommended  Values  of  Parameters  for  Tests 
(c-1)  AM  Signal  and  AM  Interferer: 


Ms  = 0.3 

•*1 = 1 

f = 0.45  W, 
bAM 


f = 0.5  W, 
AM 


- V V V £s  = fi  ■ O' 

FM  FM 


(c-2)  FM  Signal  and  FM  Interferer 


Ms  • V V "f  fs..=  *1  - 0 


AM  AM 


Og  = DJ=:  2.405 

fs  = 0.11 
*"FM 

fz  = 0.12 

FM 
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(d)  Signal  generation:  for.  some  large  number  of  positive 
values  of  J,  set 

— i 

(J)  = I ] + u ensfSrr  f ,T&  + P ] 

IG*  1 L. 


S_  (J)  = I~1  + U cos (2"  f j6  + P): 
™ S SAM  S J 


* expfj  D cos (2^  f J6)~| 

■ i?  K>  „ 


FM 


(e)  Interferer  generation:  for  some  large  number  of  positive 
values  of  J,  set 

U(J)  = |~i  + mtcos(2i*  f J6  + fl  )1 

L AM  J 

* exp|jro.lj6w  + D cos {2™  f J6)]j. 

XFM 

(f)  Noise  generation 
Calculate 

Y = "W6/2 


-r 


„ • 2 
2 sin  Y 


i~/2 


4 1/?  2 

L(l+sin  y)  ' + sin Y 


r iM^2)(i-s2)1/2]  -T 

°5  i*-2s 3 (1+s)  (i+  (i-s2) 1/2] 

Generate  a sequence  of  independent,  identically  distributed,  zero- 
mean,  unit  variance,  complex  Gaussian  variables  [g(J)}. 

Set 

n (1)  = n ( 2}  = 0. 

For  a large  number  of  positive  integers  J,  sec 

1-S  2C  SZ 

r(J)  = 2 (1-S)  n { J-l ) - — | n(J-2)  + --7-"  G ( J)  . 

1+b  1+S 
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a)  Case  selection 


The  input  data  sequences  must  be  matched  to  the  appropriate 
sequences  for  computation.  There  are  four  main  cases  which  can  be 
described  by  the  titles: 

Case  a)  Splatter  of  adjacent  channel 
Case  b)  Two-frequency  intermodulation 
Case  c)  Cross-modulation 

Case  d)  Cross-modulation  splatter  of  adjacent  channel. 

In  addition,  each  case  is  subdivided  into  two  cases  according  to 

the  assignment  of  interferer  and  noise  waveforms  to  the  internal 

data  sequences.  The  internal  data  sequences  are  labelled  Z^,Z^,Z^ 

and  z and  there  is  an  additional  labelling  variable  A which  describes 
3 

whether  the  interference  spectrum  is  centered  on  the  desired  channel 
or  the  adjacent  channel.  Table  6.3  lists  the  possible  cases. 

b)  Calculate  the  third-order  interference  product, 

Y(J>  „lj  hKj  <»<VB2'*V 

Z.tJ-S.)  32(J-!)2)  Z (J-N  ), 

and  the  linear  term 

Y, (J)  - Z ( J+K) „ 

L o 

The  required  input  information  consists  of 

(1)  the  four  data  sequences  Z , *,.,Z  rrorc  (a)  preceding, 

(2)  the  number  of  taps/ coordinate  N from  the  routine  described 
in  Section  3 

t'3)  the  coefficients  S (N, , N,, , w ) from  Section  3 

± ^ 

{4}  the  group  delay  K of  the  nonlinear  product,  from 
Section  3. 
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"“  Viw 


Note  that  the  extreme  range  of  K values  that  can  occur  is 
-129  " K < 251 

and  consequently  the  allowable  range  of  arguments  of  Y and  Yr  must 

U 

be  adjusted  accordingly. 

c)  Calculate  Y and  S as  in  Step  (f)  of  6.5.1 

d)  Calculate,  as  needed  for  output  requirements, 

(d-1)  The  total  output 

VJ)  = Yl(J)  + CtY(J). 

Required  input  parameter  is  C . 

(d-2)  The  filtered  linear  output 


Ylf(J)  = 2(M)Th(WI  - £§  \fW-2)  + ffi  Yl(J)  for  J * JQ 
with  initial  values 

Y ( J -2)  = Yt_(J  -1)  = 0, 

LF  o IiF  o 

with  J chosen  sufficiently  large  so  that  Y,  (J  ) is  defined, 
c L o 

(d-3)  The  filtered  total  output 


Y„1C.(J)  - 2(1~S)Ytf(J-1)  - Hf  Ytf(J-2)  + ffi  Yt(J)  for 


TF 
J s j 


with  initial  values 


VV2)  = yT 'Jo"  1 1 " °< 


wth  J chosen  sufficiently  large  so  that  Y,„(J  ) is  de.ined. 
o To 
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Calculate,  as  needed  for  output  requirements  the  distor- 
tion terms: 


ENV(J)  = |ytf(j)|  - |YLF(J)| 


PH  ( J)  = 


|ImlYTF/YLF}  if  'YTf'  " 2'YJ 

if  l'fTFl  * JlYLFi 


SECTION  7 


TYPICAL  OUTPUTS 

We  have  written  computer  programs  that  incorporate  all 
the  features  described  in  Sections  3,  5 and  6.  Figures  7.1  to 
7.  11  show  the  results  of  computer  runs  using  these  programs. 

The  particular  nonlinear  calculation  involved  cross- 
modulation of  a desired  signal  by  a nearby  modulated  carrier. 

We  used  t?ne  model  circuit  of  Section  5 to  create  transfer  func- 
tions and  K0.  Figures  7.1-a  and  7,.1-b  show  the  form  of  the 
computer  output  for  these  functions.  The  normalised  nonlinear 
transfer  function  is  then  computed  from  those  values  and  ie  shown 
in  Figures  7.2-a  and7.2-b.  The  normalize;:  is  also  included 

in  the  output  data  for  later  use. 

From  the  values  of  the  normalized  transfer  function  the 
optimum  model  coefficients  are  determined.  We  show  in  Figs. 

7.3-a  to  7.4-c  part  of  the  print  out  associated  vxt n the  intar- 
mediate  steps  in  this,  where  the  distortion  vs  group  delay  and 
normalized  tap  spacing  are  made  available  after  each  computation 
of  delay  line  coefficients.  The  optimum  normalized  tap  spacing 
for  this  case  proves  to  be  the  initial  value  of  D = 1/4,  and 
the  associated  group  delay  or  K is  -1.  The  computer  auto- 
matically selects  this  case  for  data  transfer  to  the  next  program 
segment,  along  with  the  correct  tap  coefficients.  This  output 
in  printed  form  is  shown  in  Fig.  7.4. 

The  generation  of  sampled  data  sequences  corresponding  to 
two  amplitude  modulated  signals  and  one  noise  waveform  is  shown 
in  Fig.  7.5. 
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The  model  parameters,  case  selection  and  J ibrary  of  three 
input  waveforms  are  then  made  available  to  the  program  segment 
that  calculates  nonlinear  outputs.  Figure  7.6  shc-v/s  the  verifies 
tion  of  the  input  selection  and  model  coefficients,  while  Fig. 

7.7  shows  the  selection  of  the  model,  inputs  according  to 


VJ) 

21(J) 

z2(j) 

z3Ut) 


“ fiIG(J) 

= SIG(J) 
- U(J) 

= U*(  T)  . 


The  nexr.  program  segment  calculates  the  third  order  inter- 
ference terms,  and  replicates  (Jl  as  the  linear  output;  these 
are  shown  in  Fig.  7.8.  The  combined  output,  with  the  weighting 

of  the  incerference  set  at  G , is  shown  in  Fig,  7.9,  This  irs — 

o 

eludes  the  proper  group  delay  N 


The  linear  and  combined  outputs  are  then  filtered  in  a 
sampled  data  second  order  Butte rworth  filter.  These  filtered 
outputs  are  shown  in  Fig,  7.10.  As  a note  of  caution,  the  first 
few  outputs  include  the  transient  response  of  the  filter,  which 
his  a response  time  roughly  equal  tc  1/D  samples.  These  first 
fevr  samples  should  not  be  used  in  subsequent  data  reduction. 


The  envelope  and  phase  distortion  are  then  calculated 
as  indications  of  the  distortion  to  be  expected  in  amplitude  or 
phase  demodulators.  These  are  shov/n  in  Fig.  7.11. 
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Computer-Generated  H and  H,  Values 
of  Data  Output  J }~ 
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Fig.  7.1b  Computer-Generated  H3  and  H.  Values?  Part  b 
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Fig.  7.8  Interference  and  Linear  Outputs 
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TinHEKTlJ  ITR-EMT  OuTPuT'*  YLE'"  ‘ 


SEQ  NO.  JO  YLF(JO) 


1 

4 

0.28400O4E 

CO 

0.0 

2 

5 

0. 5229423E 

CO 

0.  0 

3 

6 

* 0.669251 5E  ' 

CC 

"0.0 

4 

7 

0. 783141  7E 

OC 

0.0 

5 

8 

0.9174149E 

CC 

0.  C 

6 

9 

0.  107U59E 

Cl 

0.0 

7 

10 

0.  1198941E 

01 

o.c 

8 

11 

0. 1245558E 

Cl 

0.  0 

9 

12 

0.  U86191E 

01 

OYO 

10 

13 

0.  10i«446oE 

Cl 

0.  0 

11 

14 

0. 884391 5E 

CC 

0.0 

12 

15 

0.  78C3850E 

CO 

0.  c 

13 

16 

0. 7813029E 

CO 

0.0 

14 

17 

0. 8864856E 

00 

0.  0 

15 

'18 

0.  1045O37E 

oi" 

0.0 

16 

19 

J. 1132o73E 

Cl 

0.0 

17 

20 

0. 12J2067E 

Cl 

0.0 

FILTERED  TOTAL 

OUTPUT  ••  YF  T 

"SEO  NO. 

JO 

YFT(JO) 

1 

4 

0.2836743E 

00 

-0.  162384  5P-03 

2 

5 

0.  52<:544jE 

CO 

-0. 21 2519GE-03 

3 

6 

U. 66901 14E 

00 

••  0. 132523  75-33 

4 

7 

0. 7833983E 

CO 

0.701C24CE-05 

5 

8 

0.O185294E 

00 

0.3  07 4876?.-  03  " 

6 

9 

0. 1072895E 

01 

0. 6325969E-03 

7 

10 

0. 1159622E 

Cl 

0.63  J2716E-03 

8 

a 

0, 1244888E 

Cl 

0. 2^450865-03 

9 

12 

J.  1184920E 

Cl 

-0.  H35914E-03 

10 

13 

0.  1 04^>3it6E 

01 

-0.3175866^-03 

11 

14 

0. 88 J 74  78E 

00 

-0.  2*t722335-03 

12 

15 

0.  78J4j686 

00 

-0.35193885-04 

13 

16 

0.  7821930E 

CO 

0. 28J8  53  75 -03 

14 

17 

0.  68  78226E 

00 

0.57655625-03 

15 

IB 

0.  1046^16E 

01 

0.  52658°lE-03 

16 

19 

0.  lld2115E 

01 

0. 54470165-04 

17 

20 

0. 12 ->C8  76E 

Cl 

-0.3161 145E-03 

Fig.  7.10  Filtered  Outputs 
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DISTORTION  TERMS 


NO  . 

JO 

E.’iVl  JO) 

PHiJO) 

I 

4 

-0.  j320:>75E-  C3 

r 0.57176365-03 

2 

5 

-6.  39  80j  98E-  03 

-0.40639088-03 

3 

6 

-0.24014/1E-C3 

-0.22  7901  9E -03 

4 

7 

0.25&8364E-C3 

0. 8S51431E-05 

5 

8 

0.  1114488E-C2 

0.3331672E-03 

6 

9 

0.  1536369E-  C2 

0. 59C4622E-03 

7 

10 

0.68C9235E-03 

0. 52  81 924E-03 

6 

ii 

-d,<  694794E- C3‘ 

0. 1802474E-03 

9 

12 

-0.  il  C294E-  C2 

-0.  15646CCE-03 

10 

13 

-0.  . 119614E-02 

- 0. 304C65 9E-03 

11 

14 

-J.  6437302E-03 

-0.2 7554 04E-03 

12 

15 

0.  5U9o44E-C4 

-0.  5022377E-04 

13 

16 

0.8926392E-C3 

0. 3657076E-03 

u 

17 

6.  13371  71E-  C2~ 

0.  65C384C5-03 

15 

18 

0.6790161E-03 

0. 5C39885E-03 

16 

19 

-0.55789V  c -03 

0.  7587853«=-04 

17 

20 

-0.  1189232E-02 

-0.2565726E-03 

Fig.  7.11  Distortion  Terms 
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R ADC  is  the  principal  AFSC  organization  charged  with 
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reception,  ground  base d surveillance.-  ground 
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r+ocessing . 1'his  Center  provides  technical  or 
management  assistance  in  support  of  studies , analyses, 
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systems  and  related  equipment . 


Source  AFSCR  23-50,  11  May  70 


1~}% 


^ X A , ? V5A.*  * . 


